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abstr.\c:t 



Relative vorticity fields calculated from the U. S. Navy operational Global Band 
Analysis are used to relate synoptic and storm parameters to the track of tropical 
cyclones in the western North Pacific Ocean. In this preliminary study, synoptic 
patterns are developed, described and discussed from the perspective of a pattern 
recognition technique to assist the forecasters at the Joint Typhoon Warning Center, 
Guam. The focus is on track turning motions to the left and right of the persistence 
track and on tr\’ing to accurately predict the point of the turn or recurvature in relation 
to the time evolution of the vorticity patterns. The developmental sc ’e of storms 
indicates that there is potential for using synoptic patterns in the Global Band Analysis 
to guide the selection of the appropriate track aid in the 48-60 hour time range. 
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INTRODUCTION 



A. BACKGROUND 

The successful forecasting of tropical cyclone motion is one of the most difTicult 
problems in tropical meteorology. The training process for the operational forecaster 
includes a process of integrating knowledge contained in the literature and learning 
from experienced forecasters. One aspect of the training is the development of an 
ability to recognize significant patterns in the synoptic analyses, satellite imagery and 
forecast fields. The better forecasters are able to recognize the synoptic patterns that 
indicate that significant changes will occur (or not occur) in the track of the storm 
during the forecast period, or are able to recognize the synoptic patterns earlier and 
more often than the less successful forecasters. 

A significant problem at the Joint Typhoon Warning Center (JTWC), Guam, is 
the short tour length for Air Force and Naval officers assigned as forecasters. Nearly 
half of the forecaster's tour at JTWC is spent learning the job. The actual time 
required to complete the basic training phase depends on the ability and experience of 
the forecaster trainee and the number and type of storms that occur during the training 
period (McLawhorn, 1984). The short tour length at JTWC is not the only obstacle to 
meeting the COMSEVENTHFLT requirement to reduce forecast errors to 50, 100 and 
150 n mi for the 24, 48 and 72 h forecast periods respectively. The National Hurricane 
Center (NHC), .Vliami experiences forecast errors of magnitudes similar to those at 
JTWC without having to deal with the problem of forecasters subject to 15 to 24 
month rotation periods as in the military. JTWC and NHC each have several 
forecasting aids available for the development of the operational forecasts. Both centers 
are continually testing new methods developed by the research conmiunity, as well as 
making improvements to current operational models. 

It is sometimes stated that the forecaster does not need more forecast aids. 
Rather, the need is to develop the ability to better utilize the schemes that are available 
in the forecast centers (Neumann and Pelissicr, 1981; Tsui, 1984; Elsbcrry and Peak, 
1986). Each of the forecast aids has strengths and weaknesses that depend on a 
complex combination of the storm size, speed, intensity, location and the synoptic 
features in the basin. The features of importance in the western North Pacific Ocean 
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region include the subtropical ridge, long-wave troughs, short-wave troughs, the 
Northeast and Southwest monsoons, the upper-level How patterns and features, and 
otiicr tropical systems, iiach of these features is displayed to some extent in tiie 

r 

analysis and prognostic lields that are used by the forecaster, who must decide which of 
the forecast aids is the best at depicting the interaction of the environmental fields and 
the storm. In many cases, the more subtle features that are undetected by the 
forecaster under pressure to issue the forecast are later "discovered" in the post-storm 
analysis. Successful forecasts both incorporate a maximum of information into a 
warning track and help the operational customer to make prudent and timely decisions 
to save life and property at a minimal monetary cost. 

Consider the set of forecast aids available for the 12 UTC 19 August 1983 

% 

warning on Tropical Storm Dom (Fig. 1.1). This is not an atypical pattern and makes 
the selection of an accurate forecast track very dilTicult. What is required is a method 
that will incorporate the stomi and synoptic parameters into a pattern or set of 
patterns that are easily recognizable under the pressures of getting the forecast out. 
This will assist in the selection of the primary forecast aid and simplify the track 
forecasting dilemma for the operational forecaster. This feasibility study will seek the 
synoptic pattern relationships involved in the recurvature or nonrecurvature forecast. 



B. RECENT STUDIES IN SUPPORT OF THE JOINT TYPHOON WARNING 

CENTER 

One attempt at developing a technique that assists in the interpretation of the 
environmental fields is the CYCLOPS Objective Steering Model Output Statistics 
(COS.MOS) forecast aid (Allen, 1984), which uses the output from the CYCLOne 
Prediction System (CYCLOPS). The success of this technique in predicting the 
unusual track of Tropical Storm Dom is shown in Fig. I.l, COSMOS is a model 
output statistics technique applied to the CYCLOPS geostrophic steering flow aid, 
which is the most recent version of the HATIUACK/MOHATT steering flow model 
(Renard, 1968; Renard and Levings, 1969; Renard et al., 1970; Renard et al., 1972; 
Rcnard et al., 1973). The COSMOS program compares the current flow^ pattern to a 
set of seven typical patterns and calculates a forecast track that incorporates the 850, 
700 and 500 mb steering flow outputs from the CYCLOPS model. The output of a 
single track (versus separate tracks at each level) simplifies the forecaster's decision 
problem, incorporates past pattern information and allow^s the display of the 
CYCLOPS components that are used in the technique. 
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Fig. 1.1 Forecast track aids for Tropical Storm Dom 
at 12 LTC 19 Aug 19S3. 

Another recent attempt to incorporate the environmental circulation pattern into 
a forecast aid is the Colorado State University Model (CSU.M) developed bv 
-Maisumoto (19S4). This statistical technique categorizes the track forecast problem 
into one of three classes depending on the storm center location relative to liic 
subtropical ridge at 50U mb. The three classes arc defined as under the ridge, on the 
ridge and above the ridge. The predicted track is divided into three 24-hour periods 
and uses the prognostic fields as if they are analyses at the valid time. 

A series of investigations at the Naval Postgraduate School over the past three 
years has attempted to improve the use of forecast guidance at JTWC. Curry (1985) 
developed an objective technique to improve the initial location of the storm center 
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that is used for the forecast aids. Wilson (1984), Schott {1985), Williams (1986), Jones 
(1986), Meanor (1987) and W'cniger (1987) have studied various combinations of storm 
and synoptic parameters to develop information for possible ir.dusion in a deeision-tiee 
scheme for selecting the optimum track forecast guidance for a particular forecast 
scenario (Peak and Elsbcrr}', 1987). These studies have increased the knowledge base 
for the forecaster and researcher, but have not produced a straight-forward pattern that 
is easily recognizeable by the forecaster. 

Sandgathe (1987) has categorized and described the various storm-synoptic 
relationships that appear to have the greatest influence on the tropical cyclone tracks 
in the western North Pacific Ocean. He has also reviewed the 1982-1985 storm seasons 
and listed the interactions that affected the tracks of selected storms. A detailed 
description of the interactions and the forecast problems that arise due to the 
occurrence of the interaction, or a multiple interaction, are given in the report. A 
description of the interactions studied in this investigation are given in the section on 
the data set development in Chapter II. 

C. OBJECTIVES 

This study is a preliminary investig'ation of the storm and synoptic parameters 
related to tropical cyclone recurvature. The goal is to develop a way to display and 
analyze the parameters in a manner that is useful to the operational forecaster. The 
recurvature problem is chosen because of the dilTiculty experienced in accurately 
forecasting the recurvature location and time. An inconsistent performance of this 
critical forecast task often results in 72 h position errors in excess of 1000 n mi. These 
errors far exceed the COMSEVENTHFLT requirements and create major problems for 
operational commanders tasked with making decisions to implement typhoon evasion 
and preparation plans. Reduction of these track errors to levels comparable to the 
current performance for straight-moving storms would provide a credible and useful 
product to the operational decision makers. 

Several approaches to the problem have been used to discover a more optimal 
method of identifying and displaying the storm and synoptic parameters that determine 
the track of the tropical cyclone. The data fields, data retrieval and display are 
described in Chapter II, analysis techniques and results are presented in Chapter III. 
Chapter IV discusses the conclusions and possible avenues for further development and 
study. 
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II. DATA BASE DEVELOPMENT AND DISPLAY 



A. GLOBAL BAND ANALYSIS AND STORM DATA 

The wind fields in this study are taken from the United States Navy Fleet 
Numerical Oceanography Center (FNOC) Global Band Analysis (GBA). The GBA is 
the analysis used to initialize the Navy's two dynamic tropical cyclone track forecast 
models, the Nested Tropical Cyclone Model (NTCM) and the One-way Influence 
Tropical Cyclone Model (OTCM). The GBA are produced at 00 and 12 UTC each day 
on a 49 x 144 grid with grid spacing of 2.5 degrees of longitude on a Mercator 
projection. The grid extends from 40.956'"S to 59.745'N and has 360 degrees 
longitudinal coverage. The grid spacing decreases toward the poles as the meridians 
converge. At the standard latitude of the grid (22.5'^N or S), the grid distance is 257 
km. 

A complete description of the GBA can be found in the Numerical 
Environmental Products Manual (United States Naval Weather Service, 1975) and a 
review of recent developments in the Nav\‘'s tropical analysis is given in Elsberry and 
Fiorino (1985). The analysis is done at the surface, 700 mb, 400 mb, 250 mb and 200 
mb. Thermal analyses at intermediate levels provide data for the vertical coupling 
using the thermal wind relationship. The operational centers receive intermediate level 
fields calculated from analyzed levels, but these fields are not archived. When a 
tropical cyclone is present, the storm is depicted by eight wind vectors bogussed into 
the surface field around the position of the center at chart time. Data used in the GBA 
are from the surface observation network, rawindsondes, pibals, aireps and cloud 
motion vectors. If no observations are available for the current analysis, a blend of the 
12 hour old analysis with a 5% reversion to climatology is used. 

Data for tropical cyclones in the western North Pacific basin are available on an 
annual basis from the JTWC's Annual Tropical Cyclone Report (ATCR). The report 
provides a narrative description of each storm with the storm best track, satellite 
pictures, significant synoptic patterns and occasionally storm damage photographs. 
The data annex of the report gives six-hourly best track, warning and forecast 
positions; forecast error statistics; satellite, aircraft and radar fixes; best track intensity, 
aircraft wind reports and satellite intensity estimates. In this study, ATCR's for 1982, 
1983 and 1984 were used for case selection and analysis. 
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B. STORM SELECTION 



Super Typhoon Marge (18W) 1983 was chosen as the first storm to be studied 



throughout the life of the system (Fig. 2.1}. Marge was categorized by Sandgathe 
(1987. p.l9) as a step-recurver that later underwent e.\tratropical transition. A step- 
recurver is a system that is influenced by the passage of mid-latitude troughs to the 
north of the storm. The step phase of the track takes place as the storm slows and 
moves northward under the influence of a passing trough, but does not recurve. As the 
trough passes or weakens, the subtropical ridge builds westward and the storm turns to 
the left and tracks to the northwest or west under the ridge. Marge stepped twice 
before the eventual recurvature. Each time the step started, large forecast errors were 
recorded as forecasters predicted recurvature. When the storm finally did recurve, the 
forecasters unfortunately continued a track toward the northwest, which again added 
to the overall error. 



Inconsistent forecast tracks as found in the Marge case can be responsible for 
premature evacuations of aircraft and ships and the early or unnecessary’ setting of 
typhoon conditions. Both actions entail major e.xpenditures in terms of both man-hours 
and money. Table I is a summary of the average track errors for Marge compared 
with the overall errors for 1983. Although the errors and the ratio to the 1983 average 
are large, there were even larger errors on other 1983 storms in all categories except the 
72 h forecast. 

A set of storms was then compiled to test the preliminary results based on 
Marge. The goal w’as to develop a set of storms that would allow the study of varying 




TABLE 1 

TR.-\CK FORECAST ERRORS (N MI) IN 
SUPER TYPHOON MARGE COMPARED TO 1983 OVERALL ERRORS 



Marge 



1983 

Marge/ 1983 



Warning 24 h 48 h 72 h 
19 191 484 755 

16 117 259 405 

1.18 1.63 1.87 1.86 
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Fig. 2.1a Super Typhoon Marge best track (dots) and forecast tracks (stars) during 
the first step. Labels indicate date and hour (02 00 means 00 UTC 02 Nov). 
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Fig. 2.1b Super Typhoon Marge best track and forecast tracks 

during the second step. 
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Fig. 2.1c Super Typhoon Marge best track and forecast tracks 

during the recurvature. 
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storm tracks to test the pattern relationships between the synoptic Held and the storm 
envelope and to include as many storm-synoptic interactions as possible. 

Sandgacl'iC has desenbeu the track forecasting problem in the western 

North Pacilk' Ocean in terms of the storm-synoptic interactions. Mis list, which is 
ranked by need for improved forecaster understanding, includes: 



(i) Multiple cyclone interaction -- the flow created or modified by one cyclone 
influences the track of the second. 

(ii) Cyclone - midlatitude trough interaction -- the trough influences the 
cyclone track and causes a movement to the north and then recurvature if 
the trough influence is strong enough, or a step pattern if the trough 
influence is overcome by a strengthening of the subtropical ridge. Looping 
may also result from this interaction. 

(iii) Cyclone - subtropical ridge interaction — the shape and extent of the ridge 
dictates the path of the cyclone to the west or in a gradual recurvature 
around the end of the ridge. The ridge also influences the step and loop 
patterns. 

(iv) Extratropical transition -- the warm core cyclone transforms into a hybrid 
or cold-core system. 

(v) Terrain interaction -- the track and intensity of the cyclone are influenced 
by interaction of the steering flow or the storm circulation with the various 
major islands and coastal areas in the Pacific basin that have significant 
orographic features. 

(vi) Monsoon surge interaction - the tracks of some tropical storms show 
diurnal variations due to the diurnal fluctuations in the monsoon flow 
during the Southwest monsoon, which sometimes causes the storms to 
become quasi-stationary in the South China Sea and Philippine Sea. The 
Northeast monsoon also ha^ a strong influence on the track of late season 
tropical storms, and may lead to erratic tracks and large forecast errors. 

(vii) TLTT or upper low interaction — the influence of the Tropical Upper 
Tropospheric Trough (TUTT) lows on cyclone motion. 
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The 1 2 storms in the data set and the interaction categories according to Sandgathe are 
listed below (numerals following the storm name indicate the storm number and year 
and the Roman numerals correspond to tiie interaction list above): 



0) 


Bess (IIW’82) 


(ii,iii) 


(2) 


Gordon (16W/82) 


(i,ii,iii) 


(3) 


Wayne (04W/83) 


(v,vii) 


(4) 


Abby (05W/83) 


(i,ii,iii,iv,v) 


(5) 


Forrest (llW/83) 


(ii,iv,vii) 


«3) 


Ida (14W/83) 


(ii,iv,vii) 


(7) 


Lex(17W/83) 


(ii.vi) 


(S) 


Marge (18W;83) 


(iijv) 


(9) 


Ike (13 W, 84) 


(hi) 


(10) 


Thad (24W/84) 


(ii,vii) 


(11) 


Vanessa (25W/84) 


(i,ii,iv) 


(12) 


Clara (29W/84) 


(idi) 



In addition to Super Typhoon Marge, the set includes five storms that have the step 
pattern, five that recurved without first stepping and two storms that had no apparent 
interaction with the mid-latitude troughs and had straight tracks. 

C. DATA RETRIEVAL 

Wind fields are extracted from the GBA archives at 00 and 12 UTC for each of 
the storms studied. The fields include the zonal (u) and meridional (v) wind 
components at the surface, 700 mb, 400 mb and 250 mb. 

The data tapes are processed in a series of steps to allow the display of the wind 
data and derived fields in the Interactive Digital Environmental Analysis Laboratory 
(IDEA Lab) at the Naval Postgraduate School. The IDEy\ Lab consists of a 
VAX8200 mini-computer and two micro-VAX computers linked to various display 
terminals and input/ output devices. The system allows interactive display of satellite, 
observational and gridded data as well as a variety of data processing capabilities 
(Wash et aL, 1987). The first step in the process is to match the fields and storm 
tracks. A list is compiled using the ATCR best track data for each storm and is 
compared to the fields available on the data tape. A full listing of the fields available 
can be produced using the Global Band File Identification (GBFID) program found in 
Appendix A. 
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The next step is to extract the fields required for study from the Naval 
Environmental Display Network (NEDN) format tapes and create a tape in ASC H 
format for use in the IDEA Lab. Tliis is accomplished using the program GETGB.A 
listed in Appendix B. The data section of the program must be carefully prepared to 
avoid problems due to missing fields in the GBA archives. If a complete set of fields 
for a date and time are missing from the set, that date-time code entry must be omitted 
from the GETGBA data file. If only a partial set of levels or component fields is 
missing, the date-time code information should be entered, which allows the 
subroutines to extract the fields that are present. The number of records to be skipped 
to get to the desired fields must also be accurate (see documentation in Appendix B). 

The ASC II tape is used to create the data base in the IDEA Lab. The creation 
of the files for display and manipulation of the data is a two part process. Eirst. the 
tape data is read to storage on a disk using the program Global Band Tape To Disk 
(GBTTOD) found in Appendix C. The conversion of the data to a format that is 
compatible with the graphics software is accomplished using the FILE CONVERTER 
program in Appendix D. This program changes the GBA grid to a 2.5 degree latitude 
by 2.5 degree longitude grid using an interpolation subroutine, and inverts the grid 
array to make it compatible with the GEneral Meteorological (GEM) software 
(developed at NASA's Goddard Space Flight Center) used in the IDEA Lab. Once the 
data are written to a GEM FILE, the fields can be displayed and manipulated using the 
programs GEMPAK and GEMPLT. 

D. GEMPAK AND GEMPLT PROGRAMS 

The GEMPAK and GEMPLT software are used to display the data in a variety 
of ways including tabular data, grid point plots of data, contour plots, vector plots, 
streamlines and vertical profiles. The gridded data can also be manipulated vvith a 
variety of Grid Diagnostic Programs to display numerous derived fields including: 



a) 


Sums 


(ADD) 


b) 


Difference 


(SUB) 


c) 


Laplacian 


(LAP) 


cl) 


Layer Average 


(LAV) 


c) 


Time Derivative 


(DDT) 


0 


Time Difference 


(TDF) 


g) 


Time Average 


(TAV) 


h) 


Divergence 


(DIV) 


i) 


Relative Vorticity 


(VOR) 
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j) Absolute Vorticity (AVOR) 

k) Gradients (GRAD) 

and a number of other functions found in the GUMPAK L’ser's Guide i desJardin ci a!.. 
19S6) Some of the diagnostic grids can be created "on the lly" and be directly displayed 
or sent to a hard copy output device. More complicated functions require that 
intermediate grids be created before further processing. In either case, an option exists 
to save the diagnostic grids in the same GEM file as the original storm data in date- 
time order. 

The display terminals give multi-color depictions of the fields including 
geography of various resolutions and projections. Fields can be overlayed on the 
screen for analysis and can be reproduced in hard copy on laser printers or full color, 
on the VERSATEC plotter. 
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III. DATA ANALYSIS 
A. ANALYSIS PARAMETERS 

Analysis of the fields is based on significant storm and synoptic parameters and 
the changes with time in the synoptic fields. The storm parameters used in the study 
are those used by the operational forecaster to request and initialize the forecast aids 
generation programs at FNOC. The parameters used in the study (Table 2) differ 
slightly from the combined aids request (CARQ) as the storm position information is 
reduced from 48 to 36 h and the intensity and size data are increased to 36 h of history. 
The storm size information is taken from the waniins mcssasies usina the wind radii for 

W W ^1.. 

30, 50 and 100 kt. The values relating to the outer closed isobar as required in the 
forecast aids request were unavailable. 



table 2 

STORM PAR/\METERS 

1. Month of year, day of month 

2. Position 

a. Latitude, longitude 

b. Fix platform (aircraft, satellite, radar, synoptic or extrapolation) 

c. Position and direction change - movement past 12, 24, 36 h 

3. Intensity 

a. Surface wind (kt), MSLP (mb) 

b. 700 mb winds (kt), height (m) 

c. Change in above parameters in the last 12, 24, 36 Ii 

4. Storm size 

a. Radii of 100, 50 and 30 kt winds 

b. Change in radii in the last 12, 24, 36 h 
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The synoptic parameters (Table 3) include those required for the CARQ as well 
as features used in other forecast techniiiues or deemed significant through 
developmental research Iindings. Tiie position and intensity oi' the subtropieai ridge 
and mid-latitude trough are included as they are part of the CARQ values, and the 
history change values are extended to 36 hours. The polar and subtropical jets are 
considered in light of the Typhoon Acceleration Prediction Technique (TAPT) 
developed by Weir (1982). The positions of the short wave trough north of the 
tropical cyclone were added to the list after the initial analysis had begun, as small 
deviations in the track were often found to coincide with the short-wave trough 
passage. As in the storm parameters, the changes in the synoptic parameters are 
recorded for the past 36 h. 

B. STREAMLINE FIELD ANALYSIS 

The initial analysis of the data involved a detailed study of Super Typhoon 
Marge using the streamline version of the GBA wind fields (as used at JTWC) and the 
National Meteorological Center (NMC) Tropical-Mercator surface and upper-air 
analyses. This approach proved difficult for several reasons. The NMC analyses were 
not available for all time periods, and the mid-tropospheric chart was for the 500 mb 
vice the 400 mb level in the GBA series. This led to problems in developing a 
continuous history and determining accurate height center locations corresponding to 
the wind field features. The GBa\ streamline output presented inconsistencies in the 
depiction of some small centers, and often showed cyclones and anticyclones with 
reverse flow patterns. While the bearings from the storm center to the synoptic 
features could be determined in a consistent manner, the range determinations were 
often suspect due to scale changes across the Mercator grid. After completing the 
analysis of Marge using this approach, several conclusions were made: 

a. The process requires too much time to be useful in a forecast situation; 

b. The results arc not clear due to the large number of variables involved; and 

c. Inconsistencies from chart to chart reduce the value of the information. 

C. VORTICITY FIELD ANALYSIS 
1. Vorticity Field Presentations 

The storm and synoptic parameters are important and require a presentation 
format that will incorporate the variations in the parameters into a more easily 
interpreted form. This idea led to the use of the vorticity fields (Fig. 3.1), which are 



15 



'TABLi: : 

S \'N0 P'l' 1 C PA FL\ M liin RS 



1. Subtropical ridge (Surface, 700 and 400 mb) 

a. range and bearing to the maximum height of the ridge 

b. range and bearing to the lowest height or break 

c. height magnitude range and extent of the ridge 

2. Long-wave trough (Surface, 700 and 400 mb) 

a. ranae and bearing to the trough axis at 35 deg N to the west 

WWW w 

b. height of the trough 

c. closing rate of trough to the storm center 

3. Short-wave trough 

a. range and bearing to trough to the NW-NE of the storm center 

b. closing or opening rate 

4. Upper-level trough 

a. TUTT location 

b. closed cell locations and intensity 

5. Subtropical and polar jets 

a. range and bearing to the jet max XW-NE of the storm center 

b. TAPT pattern type (Weir, 1982) 

6. Second, third or fourth storm in the basin 

a. bearing and range to each storm 

b. closing rates 

c. pressure and wind speed differentials 

7. Topographic influence 

a. range and bearing to the nearest influential land feature 

b. flow direction and speed over the feature 

c. closing or opening rate between the storm and feature 



presently not available to the JTWC forecaster. The bogussed storm and the analyzed 
features in the GBA wind field are displayed well in the relative and absolute vorticity 
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(lelds using die diagnostic capabilities of' the GEMPAK software in the IDEA Lab. 
The conversion to vorticity eliminated the reverse How around small features that 
appeared in the streamlines and ailowed better analysis of tlie details in the major 
features. The subtropical ridge in the vorticity fields gives details of the relative 
strength of dill'erent areas of the ridge that are not represented well in the streamline 
chart. The intensity of the synoptic troughs is also better shown in the contoured 
vorticity chart. The packing of the contours around the low centers shows the relative 
strengths of the centers as well as the changes in intensity when viewed in a time 
sequence. Thus, the vorticity presentation conveys important information to the 
forecaster at a glance without requiring lengthy data tabulation. The smaller scale 
storm and synoptic features appear to be integrated into the vorticity fields and 
presented through the variations found in the contour patterns between the major ridge 
and trough depictions. 

It is hypothesized that the presence and location of a weakness in the ridge is 
very important when trying to predict the longitude of a turning or recurvature event. 
Perhaps an improvement in the prediction of the recurvature point could result from 
using the vorticity charts in conjunction with the objective track forecast aids. 

DeMaria (19S5) suggests that the storm motion can be related to the absolute 
vorticity gradient of the steering current with motion components in the direction of 
the gradient and cross-gradient. Thus the analysis of the vorticity fields was done 
initially using the gradients of relative and absolute vorticity. The model results of 
DeMaria were not evident in the operational fields, perhaps because of errors and 
inconsistencies in the operationally-analyzed fields that are then magnified by taking 
gradients. Since the error characteristics of these vorticity fields are not well known, 
the forecaster must determine whether a vorticity feature represents a real circulation 
or is caused by data inconsistencies. 

The relative and absolute vorticity fields contain several dominant patterns 
that are found to vary gradually in time as the interactions among the tropical storm 
and the adjacent synoptic features strengthen and weaken with time. The relationship 
of the cyclone track to. the synoptic patterns revealed in this study (Fig. 3.2) are: 

(i) CHANNEL — linear pattern with relatively straight relative vorticity 
contours on each side of the storm that results in a "cork in a stream" type 
pattern (Fig. 3.2a). That is, the past and future storm positions tend to lie 
within the same "channel" of relative vorticity. 
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Fig. 3.1a Streamline chart at 400 mb at 12 UTC 31 Oct 1983. 
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Fig. 3.1b Relative vorticity (xlO'%*^) chart at 400 mb at 12 L’TC 31 Oct 19S3. 
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ui) NODE - bulges in Uie contour(s) around the storm center that extend 
OLit'.vard from the average circumference of the contour around the center, 
or extend northward from the average latitude of a zonally-oriented 
contour to the north of the center (Fig. 3.2b). Primary’ nodes are defined to 
be along the major axis of an elliptical vorticity contour that encloses the 
storm, or along the persistence track. That is, the tropical cyclone appears 
to track toward a primary node. However, secondary nodes are also 
monitored because they may become the major axis of the elliptical center, 
and thus be related to the subsequent storm direction changes with time. 

(iii) CLOSED CONTOURS -- vorticity (positive and negative) centers adjacent 
to the storm center that potentially may influence the long-term track (Fig. 
3.2c). 

(iv) PINCHED or BROKEN RIDGE — areas of contour deformation in the 
ridge north of the center that indicate weaknesses in the ridge (Fig. 3. 2d). 

(v) STORM-STOR.VI INTERACTION - pattern transition from independent 
closed contours enclosing separate storm systems to a continuous vorticity 
contour enclosing two or more storm centers (Fig. 3.2e). 

(vi) RESULTANT VECTOR PATTERN -- a closed vorticity contour with a 
primary and a secondary’ node with axes approximately 80-120 deg apart 
(Fig. 3.2f)- In some cases, a resultant vector from these two nodes can be 
constructed that appears to have very/ good correlation with the track at 
60-72 h. 

2. Proposed Synoptic Analysis Technique 

The proposed analysis technique involves a recognition of the dominant 
synoptic pattern from the above list. Changes in the synoptic pattern with time are 
then recorded and are related to expected changes in the storm track. The orientation 
and deformations in the zero contour of relative vorticity is the key to the analysis, as 
it is common to all plots and usually encloses the storm center. Based on George 
(1975), Chan and Gray (I9S2). Chan (1984) and Chan (1985), vorticity fields at 700 
mb, 400 mb and the layer average from 700-400 mb are examined to define synoptic 
pattern relationships to the track of the storm. The layer average is more consistent in 
time than the 700 and 400 mb fields alone. The Utver averaae dves a smoother contour 
field that seems to change more conservatively through time than either of the 
individual levels, as suggested by Chan and Gray (1982). The relative vorticity fields 
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Fig. 3.2a Relative vorticity contour - channel pattern. 




Fig. 3.2b Relative vorticity contour - node pattern. 




Fig. 3.2c Closed relative voriicity contour pattern. 




Fig. 3. 2d Pinched or broken ridge contour pattern. 




Fig. 3.2e Storm - storm interaction contour pattern. 




Fig. 3.2f Dual node-derived resultant vector pattern. 
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are chosen because the patterns of interaction among the storm and synoptic features 
are clearer with the earth's vorticity removed. 

In the development phase of the research, a set of relative vorticity charts was 
plotted for the duration of Super Typhoon Marge. The best track for the 12 h old 
through the 72 h future positions were plotted on the charts to allow comparison of 
the vorticity pattern and the storm track. Unfortunately, the one chart set that was 
not available (12 UTC 05 Xov 1983) in the Marge series was 6 h after the time of 
recurvature. Table 4 eives the svmbols used on the storm charts in the remainder of the 
study. 



TABLE 4 

STOR.M CHART LEGEND 

Best track positions (12 h intervals) 
Current position # 

12 h old position ■ 

Future positions • 

Forecast track positions (12 h intervals) # 

Positive vorticity contours 

Negative vorticity contours 



Each chart is analyzed by first checking around the storm center location to 
determine the dominant contour pattern. In Fig. 3.3a, the first reported position of 
Typhoon Marge is located on the southern edge of a closed zero contour at 400 mb. 
The major pattern in the contours is a primary node in the zero and -2 contours, which 
has an orientation of approximately 310 degrees. Secondary nodes are found at 340 
deg in the 0 contour and 325 deg in the -2 contour. Three distinct centers are found in 
the subtropical ridge, which is well to the north of the storm at this time. Notice that 
lines from the storm center through the primary^ and secondary nodes appear to pass 



24 



through the breaks between the centers in the subtropical ridge. The vorticity contours 
on the northern side of the ridge appear to have little relationship to the tropical 
cyclone as they are influenced more by the positive vorticity centers associated with 
two extratropical cyclones to the north of 40®N. 



BEST TRACK 


TABLE 5 

- RELATIVE VORTICITY PATTERN 


ANGULAR CORRELATION EVALUATION CRITERIA 


Included Angle (a) 


Quality Assigned 


Points 


0 < a < 5 


Excellent 


6 


5 < a < 15 


Good 


5 


15 < a < 30 


Fair 


4 


30 < a < 60 


Poor 


2 


60 < a 


None 


0 



To describe the correlation between the contour pattern and the best track, a 
comparison of the angle between the contour feature and the best track position is 
made. A point value is assigned according to the magnitude of the angle between the 
pattern axis and the best track. Table 5 gives the breakdown of points and a qualitative 
description for the ranges of the included angle (a) considered. 

In applying the evaluation scheme to Fig. 3.3a, the closed zero contour that 
the center lies on is not considered. The orientation of the primary node in the open 
zero contour is evaluated as excellent (6 points) at 42 h and the secondary node is 
judged to be fair (4 points) at 42 h. These times are estimated by drawing an arc from 
the storm position and the maximum extent of the contour pattern to the best track, 
and interpolating along the best track at the intersection point (if the time is greater 
than 72 h, a value of 72 is assigned). For example, the primary and secondary nodes on 
the -2 contour are evaluated as excellent and good respectively at 72 h. 

The vorticity contours in the 700 mb and 400:700 mb layer average fields (Fig. 
3.3b and c) have less detail near the storm, although both have a node to the 
northwest in the zero contours. The 700 mb node in the zero contour is rated excellent 
at 58 h. The small +2 closed contour is evaluated as excellent at 12 h. The layer 
average (LAV) node in the zero contour is rated to be excellent at 48 h. Notice that a 
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Fig. 3.3a Relative vorticity (xlO^s"^) at 400 mb at 12 LTC 30 Oct 1983. 
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Fig. 3.3b Relative vorticity (xlO^s*^) at 700 mb at 12 L’TC 30 Oct 1983. 




Fig. 3.3c Relative vonicity (xlO-^s‘*) in LAV at 12 UTC 30 Oct 1983. 
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distinct break in the ridge is defined by the -2 contour north of the 72 h best track 
position. 

The movement and development of these basic patterns is followed through 
the next eight days in 12 h intervals to illustrate the progression of these synoptic 
pattern-track relationships. The objective is to determine if the synoptic pattern 
changes precede storm track changes consistently, which would allow them to be useful 
for forecasting the track change. 

As development occurs in the next 12 h, a new +2 contour is found around 
the storm (Fig. 3.4a). Although the orientation of the zero contour nodes change little 
in 12 hours, the -2 node has moved slightly south and closer to the storm center in 
response to a strengthening of the subtropical ridge to the northwest. Notice that the 
change in the -2 node correlates with the direction change in the best track at 60-72 h. 
The strenghtening of the ridge to the northeast seems to be linked to the passage of a 
vorticity center (45®N, 170®E) moving to the east-northeast, with strengthening of the 
ridge following passage of the +12 center and weakening in proximity of the +6 
center northwest of Marge. The 700 mb field (Fig. 3.4b) has changed little, other than 
a widening of the ridge in the meridional direction and the enclosure of the storm 
center by a +2 contour. Similarly, the LAV (Fig. 3.4c) has an enclosed vorticity center. 
Although the ridge has strengthened, the break in the ridge near 25°N, 140°E has 
closed. This again suggests the storm may not continue the northerly track. 

The primary (secondary) node on the zero contour at 400 mb becomes more 
(less) pronounced in the next 12 h (Fig. 3.5a). However, a large secondary node 
develops on the -2 contour. Although the new node on the -2 contour is larger than 
the one that has been tracked previously, it is considered as secondary due to the lack 
of past history for the feature. Furthermore, the new node does not correlate with the 
primary zero contour node or the storm track history well. At 700 mb (Fig. 3.5b), an 
expansion of the vorticity center occurs around the tropical cyclone, a closed zero 
contour develops and some weakening of the ridge occurs. A distinct change has taken 
place in the LAV chart (Fig. 3.5c) as a primary node in the zero contour has developed 
to the northwest of the storm center. This node is rated as a good correlation (5 
points) at 60 h. However, the subtropical ridge has also widened slightly and a break 
has redeveloped to the north of the storm center in response to the mid-latitude trough 
pattern. The forecaster would have to evaluate whether this break in the ridge would 
signal continued recurvature, or whether the development of a strong primary node to 
the northwest would indicate a turn to the west. 
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Fig. 3.4a Relative vorticity (xlO^s'*) at 400 mb at 00 LTC 31 Oct 1983. 
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Fig. 3.4b Relative voriicity (xIO^s*^) at 700 mb at 00 UTC 31 Oct 1983. 
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Fig. 3.4c Relative vorticity (xlO^s"^) in LAV at 00 UTC 31 Oct 1983. 
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Fig. 3.5a Relative voriicity (xlO-^s'h at 400 mb at 12 UTC 31 Oct 1983. 
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Fig. 3.5b Relative vonicity (xlO'^s*^) at 700 mb at 12 UTC 31 Oct 1983. 
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Fig. 3.5c Relative voniciiy (xlO^s**) in LAV at 12 UTC 31 Oct 1983. 
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At 00 UTC 01 Nov 83, the 400 mb relative vorticity field (Fig. 3.6a) has slight 
changes including a +2 closed contour around the tropical storm center and 
strengthening of the ridge to the north of the center. Although continued gro\nh of the 
secondary node in the -2 contour is observed, the tropical cyclone appears to be 
embedded in the primary node that in this case is associated with sustained motion 
toward the northwest. In Fig. 3.6b, the continued expansion of the vorticity center of 
the storm is evident as well as the appearance of a pronounced node in the zero 
contour to the northwest. The large vorticity gradient oriented NW to SE to the 
northwest of the storm would seem to suggest continued movement to the northwest. 
The LAV pattern (Fig. 3.6c) also shows slight changes. Although the subtropical ridge 
break has closed, a definite pinched ridge pattern remains, and the meridional 
separation of the -2 contours is reduced. The zero contour has changed little and a + 2 
contour has enclosed the center. The -2 contour changes would seem to indicate that 
recurvature should be anticipated at this time. However, the zero contour at all levels 
continues to display a major axis in the northwest-southeast direction, which would 
indicate movement to the northwest. 

At 12 UTC 1 Nov 83, a significant change occurs at 400 mb (Fig. 3.7a) as the 
zero contour now has a pronounced secondary node similar to the node in the -2 
contour. Thus, the direction of the major axis of the pattern is less clear than before. 
However, these nodes are quite far north of the tropical cyclone. Also, note that the 
major axis of the + 2 closed contour around the center remains in the northwest- 
southeast orientation. Based on the persistence of the primary node and the recent 
track of the storm, no change is made in the primar>^ node designation at this time. 
The 700 mb chart (Fig. 3.7b) shows a direct link to the higher latitudes has developed 
in the form of a narrow channel through the ridge near 22®N, HO^E. The LAV chart 
(Fig.3.7c) has a more conservative change in the contour pattern with only a gradual 
growth of the secondary node on the zero contour. Although the ridge to the north of 
the center is quite b^^oad, the ridge to the west-northwest appears to be narrowing in 
response to a positive vorticity center now over Honshu and moving to the east. This 
center eventually will create the path to the mid-latitudes that STY Marge will follow 
in recurvature. In hindsight, the key to the present situation is to recognize that the 
opportunity for recurvature into the weakness in the ridge will not occur. Only a feint 
to the north actually occurred, presumably because the primary vorticity channel 
oriented NW-SE was a dominant factor rather than the possible influence of the more 
distant secondary 
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Fig. 3.6a Relative vorticity (.xlO^s**) at 400 mb at 00 UTC 01 Nov 1983. 
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Fig. 3.6b Relative vorticiiy (xlO-%'*) at 700 mb at 00 UTC 01 Nov 1983. 
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Fig. 3.6c Relative voriicity (xlO^s*^) in LAV at 00 UTC 01 Nov 1983. 
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Fig. 3.7a Relative vorticity (xlO^s"^) at 400 mb at 12 UTC 01 Nov 1983. 
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Fig. 3.7b Relative vorticity (xlO“s'^) at 700 mb at 12 LTC 01 Nov 1983. 
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Fig. 3.7c Relative vorticity (xlO^s*^) in LAV at 12 UTC 01 Nov 19S3. 
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The next chart sequence (00 UTC 02 Nov 83) shows a return to a definite 
northwest track orientation in the 400 mb and LAV charts. The sharp secondary node 
is no longer present in the zero contour at 400 mb (Fig. 3.8a). The 700 mb chart (Fig. 
3.8b) shows the further development of the channel pattern near 20®N, HO’^E, which 
connects with the +2 vorticity center south of Honshu. Although a clear opening in 
the ridge is evident to the northwest in the LAV (Fig. 3.8c), the small positive center 
detected on the previous LAV chart has moved to the south of Honshu. Nevertheless, 
a weak cell in the subtropical ridge remains directly between the mid-latitude trough 
and the tropical cyclone. 

At 12 UTC 02 Nov 83, strong indications of recurvature are seen in all three 
charts (Fig. 3.9). In Fig. 3.9a (400 mb), the primary node on the zero contour is clearly 
to the west-northvvest. This primary node on the zero contour shows excellent 
correlation to the best track (6 points) at 36 h. However, a definite shift is seen in the 
primary^ node on the -2 contour in response to the continued strengthening of the 
positive vorticity center near 32’^N, 142‘*E (good rating 5 points at 72 h). This change 
in the -2 contour nearly breaks the ridge at 24°N, 135’"E. By contrast, the secondary 
node on the -2 contour shows little change at this time. Recurvature is expected late in 
the period, which is confirmed by the good (5 points) rating at 72 h. The cyclone center 
is in a channel of positive vorticity at 700 mb (Fig. 3.9b), with the channel wall (zero 
contour) to the north nearly describing the eventual recurvature path 66 h before the 
recurvature takes place. The LAV (Fig. 3.9c) indicates the impending recurvature in 
another way. With the eastward retreat of the -2 contour that previously separated the 
positive vorticity center to the north from the tropical cyclone center, the distance 
between the zero contours separating the centers is reduced as the positive center south 
of Honshu moves to the east-northeast. Although a clear pathway to the mid-latitudes 
appeared to develop in this case, the early portion of the track was toward the west- 
northwest along the primary node in the zero contour. 

The changes in the individual levels (Fig. 3.10a and b) at 00 UTC 03 Nov 83 
are somewhat confusing. At 400 mb (Fig. 3.10a), the ridge strengthens substantially, 
which leaves only a narrow break to the northeast of the storm center that is aligned 
with the secondary node in the zero contour. A reduction of the primary node in the 
zero contour is also evident. The channel at 700 mb (Fig. 3.10b) that had persisted for 
36 h has closed, which leaves a weakness in the ridge similar to that at 400 mb. The 
LAV chart (Fig. 3.10c) reflects the alignment of weaknesses at the upper and lower 
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Fig. 3.8a Relative vorticity (xIO^s'^) at 400 mb at 00 UTC 02 Nov 1983. 
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Fig. 3.8b Relative vorticity (xlO“s‘*) at 700 mb at 00 UTC 02 Nov 1983. 
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Fig. 3.8c Relative vorticiiy (xlO^s*^) in LAV at 00 UTC 02 Nov 198 
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Fig. 3.9a Relative vorticity (xlO^s’^) at 400 mb at 12 UTC 02 Nov 1983. 
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Fig. 3.9b Relative vorticity (xlO^s’^) at 700 mb at 12 UTC 02 Nov 1983. 
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Fig. 3.9c Relative vorticity in LAV at 12 UTC 02 Nov 19S3. 
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bounds of the layer, and also indicates a flattening of the primary node on the 0 
contour, as well as an increase in the secondary* node amplitude. 

All three charts at 12 UTC 03 Nov S3 show strong recurvature tendencies as 
the indications of 24 h previous are reconfirmed. The 400 mb and LAV charts (Fig. 
3.11 a and c) both show substantial increases in the node to the northeast of the 
cyclone center. At 400 mb, the primary node still shows excellent correlation with the 
best track (6 points) but the time factor has been reduced to about 30 h. Not only have 
the secondary nodes in the 0 and -2 contours increased in amplitude, they are aligned 
with each other along an axis that is nearly parallel with the post-recurvature track of 
STY Marge. The LAV chart (Fig. 3.11c) shows a greater flattening of the primary node 
in the zero contour (rated as excellent or 6 points at 36 h) and a marked increase in the 
amplitude of the secondary node. A line drawn from the cyclone center through the 
secondary node and the greatest southward extent of the zero contour to the north of 
the ridge also correlates well with the post-recurvature storm track. At 700 mb (Fig. 
3.11b), the return of the southwest-northeast channel defined by the zero contours is 
clearly indicated. The axis of the channel correlates well with the patterns in the other 
two layers and the post-recurvature track. Consequently, the key question at this stage 
is not whether recurvature will occur, but when? At the present speed and direction of 
the storm, it would arrive at the flattened portion of the nearest zero contour in about 
30 h. A turn toward the north should be expected at that time. 

At 00 UTC 4 Nov 83, the designation of a new primary node must be made 
on all three charts. Although the present track of the cyclone is still aligned with the 
old primary node, the positive vorticity area around the cyclone center has a distinct 
shift in the orientation of the major axis toward the northeast. Thus, the old secondary 
node becomes the new primary node. Although the 400 mb and LAV charts (Fig. 3.12a 
and c) show similar patterns with large nodes to the northeast, an additional pattern is 
also present. The -4 contour on the 400 mb and the -2 contour on the LAV display a 
pinched ridge pattern (Fig. 3. 2d) that aligns well with the post-recurvature path. The 
700 mb chart (Fig. 3.12b) includes an extended channel pattern far to the northeast 
that encompasses the "wake" of the positive vorticity center that started as a small 
center south of Korea (Fig. 3.6c) and tracked to the east-northeast to the north of STY 
Marge. The low' center passed to the north of STY Marge at 12 UTC 02 Nov 83 (Fig. 
3.9c) and passed from the chart area after 12 UTC 03 Nov 83. That is, a frontal trough 
remains at 700 mb along the expected recurvature path. This feature splits the 
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Fig. 3.10a Relative vorticity (xlO^s**) at 400 mb at 00 UTC 03 Nov 1983. 
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Fig. 3.10b Relative vorticity (xl0~s*^) at 700 mb at 00 UTC 03 Nov 1983. 
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Fig. 3.10c Relative vorticity (xlO^s'^ in LAV at 00 UTC 03 Nov 1983. 
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Fig. 3.1 la Relative vorticity (xlO^s**) at 400 mb at 12 LTC 03 Nov 1983. 
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Fig. 3.11b Relative vorticiiy (xlO~s‘^) at 700 mb at 12 UTC 03 Nov 1983. 



55 




Fig. 3.11c Relative vorticiiy (xlO^s'^) in LAV at 12 LTC 03 Nov 1983. 
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subtropical ridge around 25°N and the mid-latitude ridge that extends southward from 
40*N, 160°E. 

By 12 UTC 04 Nov 83, an extensive node of positive vorticity has penetrated 
into the subtropical ridge during the past 24-36 h. In Fig. 3.13a (400 mb), the node axis 
continues to be nearly parallel to the best track after recurvature with a slight right 
bias. Thus, the recurvature track appears to punch directly through a rather substantial 
ridge at 400 mb. However, Marge is clearly approaching the ridge axis at 700 mb (Fig. 
3.13b). The zero contour around Marge at 700 mb and the one associated with the 
mid-latitude positive vorticity band are much closer than those at 400 mb or previously 
at 700 mb (Fig. 3.12b). This indicates that the subtropical ridge is narrowing at the 
lower levels more than at 400 mb at this time. Both lower and upper-level charts need 
to be examined in these baroclinic (vertically-sheared) conditions. The LAV chart (Fig. 
3.13c) has a combination of a broad node to the northwest and a narrowing of the 
subtropical ridge. The LAV node axis is not as clearly correlated with the post- 
recurvature track as is the 400 mb depiction. This poor correlation may be due to the 
layer averaging routine that does not use a mass weighted scheme, but rather a simple 
averaging of the two input level values. That is, the 700 mb level should have much 
more weight in the LAV fields because it represents a deeper layer of the atmosphere 
that is steering the storm. 

Over the next 12 h, the trends continue in all three charts (Fig. 3.14 a, b and 
c). A continued reduction of the distance separating the zero contours of the ridge is 
seen. The next 12 h period is not available as the 12 UTC 05 Nov 83 fields are missing 
from the archived GBA data. Thus, the transition at the time of recurvature between 
00 and 12 UTC cannot be observed for Marge. 
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Fig. 3.12a Relative vorticity (xlO^s'^) at 400 mb at 00 UTC 04 Nov 1983. 
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Fig. 3.12b Relative vorticity (xlO^s’*) at 700 mb at 00 LTC 04 Nov 1983. 
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Fig. 3.12c Relative vorticity (xlO^s’^) in LAV at 00 UTC 04 Nov 1983. 
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Fig. 3.13a Relative vorticity (xlO^s"^) at 400 mb at 12 UTC 04 Nov 1983. 



61 



I 




Fig. 3.13b Relative vorticity (xlO^s*^) at 700 mb at 12 UTC 04 Nov 1983. 
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Fig. 3.13c Relative vorticity (xlO^s**) in LAV at 12 UTC 04 Nov 1983. 
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Fig. 3.14a Relative vorticity (xIO'^s"^) at 400 mb at 00 UTC 05 Nov 1983. 
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Fig. 3.14b Relative vorticity (xlO^s*^) at 700 mb at 00 UTC 05 Nov 1983. 
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Fig. 3.14c Relative vorticity (xlO^s'b in LAV at 00 UTC 05 Nov 1983. 
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A major development of a positive vorticity center at 400 mb (Fig. 3.15a) 
occurs over southern Japan at 00 LTC 06 Nov 83. Consequently, only a very narrow 
ridge separates the tropical cyclone and the mid-latitude circulation. At 700 mb (Fig. 
3.15b), a +2 contour surrounds both the tropical and mid-latitude systems with a large 
node in the zero contour on the east side of a north-south channel. The +2 closed 
center in the node has excellent correlation with the cyclone's final best track position. 
The LAV (Fig. 3.15c) includes two distinct centers in a positive vorticity channel with 
about 300 n mi separation of the + 2 contours around the centers. A reflection of the 
700 mb node is present in the -2 contour to the east. Again, a more representative 
vertical weighting in the LAV would probably display the recurvature conditions more 
accurately. 

Over the next 24 h, the track of STY Marge appears to be determined by that 
of the mid-latitude system tracking to the east-northeast. At 400 mb (Fig. 3.16a and 
3.17a), the tropical center seems to be moving slightly faster than the mid-latitude 
system. In the 700 mb charts (Fig. 3.16b and 3.17b), only one vorticity center is 
depicted. Notice the extreme vorticity gradient between the tropical cyclone and the 
ridge to the east. The best track tropical cyclone position is to the southeast of the 
large center. The LAV charts show a blend of the two levels as expected (Fig. 3.16c 
and 3.17c). 
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Fig. 3.15a Relative vorticity (xlO^s’M at 400 mb at 00 UTC 06 Nov 1983. 
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Fig. 3.15b Relative vorticity (xlO^s*^) at 700 mb at 00 UTC 06 Nov 1983. 
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Fig. 3.15c Relative vorticity (xlO^s"^) in LAV at 00 UTC 06 Nov 1983. 
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Fig. 3.16a Relative vorticity (xlO^s’^) at 400 mb at 12 UTC 06 Nov 19S3. 
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Fig. 3.16b Relative vorticity (xlO^s*^) at 700 mb at 12 UTC 06 Nov 1983. 
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Fig. 3.16c Relative vorticiiy (xlO^s’*) in LAV at 12 LTC 06 Nov 1983. 
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Fig. 3.17a Relative vorticity (xlO^s'^) at 400 mb at 00 UTC 07 Nov 1983. 
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Fig. 3.17b Relative vorticity (xlO^s'b at 700 mb at 00 UTC 07 Nov 1983. 
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Fig. 3.17c Relative vorticity (xlO'^s*^) in LAV at 00 UTC 07 Nov 1983. 
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3. Summary of the Relative Vorticity Analysis of STY Marge's Track 

The development of the node pattern in the relative vorticity contours and the 
change in the pattern with time appear to give reasonable indications of the three 
critical changes in the track of STY Marge. The first step in the track occurred at 12 
UTC 31 Oct 83 and was the smallest of the track changes. Each of the 400, 700 and 
LAV chart series indicated that the subtropical ridge was strong. The primary nodes in 
the vorticity patterns showing excellent correlation 18-30 h beyond the 20 deg left turn 
in the track. The strongest indicator was the 400 mb pattern, which provided the best 
indication of the conditions during the first step in the storm track. 

The weakness in the ridge that accompanied the bending of the storm track to 
the risht after 12 UTC on the 31st is evident in the 400 mb and LAV charts, but the 
pattern does not appear clearly at 700 mb. The ridge at 700 mb narrows to some 
degree but remains intact to the north of the cyclone center. The ridge break that 
develops to the northeast of Luzon is well to the west-northwest of the center. In this 
case, the pattern at 700 mb would indicate west-northwest movement rather than 
recurvature, about 48 h prior to the left turn at 12 UTC on 2 Nov 83. The 400 mb 
chart 24 h prior to the step in the storm track (Fig. 3.7a) could lead to a serious track 
forecast error if it was used without the other two charts. The sharp north-northeast 
node would seem to indicate that recurvature should take place at about 140°E. 
However, the 700 mb chart indicates a relatively strong ridge at MO^E with a weakness 
in the ridge to the west near HO^E. The average of the ridge weaknesses at the two 
levels would be at 135®E, which coincidentally would have been an excellent forecast as 
STY Marge eventually did recurve at about 134*E. 

The critical recurvature period appears to be well represented in the relative 
vorticity fields as the 400 mb node pattern and the 700 mb channel pattern began to 
align in the vertical after the second track step. The transition of the primary- node 
pattern from one oriented toward the NW to become oriented with the previous 
secondary node started 72 h prior to the recurvature, and was nearly complete 36-48 h 
prior to the turn to the northeast. This would seem to provide a useful forecast aid if 
such patterns appeared consistently. Finally, orientation of the node and channel axes 
gave an indication of the post-recurvature track. The question is whether the patterns 
observed in the relative vorticity fields of STY Marge are evident in other tropical 
cyclones. Table 6 gives the dominant vorticity pattern, rating value and time factor for 
each 12 h period for STY Marge. 



77 




78 




D. APPLICATION OF THE ANALYSIS TECHNIQUE TO OTHER STORMS 
1. Dominant Vorticity Patterns 

The relative vorticity contour pattern analysis described above is applied to 
the other 11 storms in the data set. It is not possible to describe these cases in the 
detail given above for STY" Marge. Dominant and secondary contour patterns (Fig. 
3.2) are identified and the angular difference (Table 5) between the pattern and the best 
track is measured, and the time interval over which the technique would provide useful 
guidance is estimated as above. Table 7 is a summary of the relative vorticity contour 
patterns observed in the storms. 



TABLE 7 

RELATIVE VORTICITY CONTOUR PATTERNS OBSERVED 



Storm 


Dominant Patterns 


Secondary Patterns 


Marge 


Node 


Ridge deformation, 
Channel 


Bess 


Node, 

Channel 


Ridge deformation. 
Resultant vector 


Gordon 


Node, 

Channel 


Ridge deformation, 
Storm- storm (Faye) 


Forrest 


Ridse deformation. 
Node 


Closed centers. 
Resultant vector 


Lex 


Node 


Channel 


Ida 


Node 


Closed centers. 
Resultant vector 


Clara 


Ridge deformation. 
Node 


Closed centers. 
Storm-storm (Bill) 


Thad 


Ridge deformation, 
Channel 


Closed centers 


Vanessa 


Node, 

Channel 


Storm-storm (Thad, Warren) 


Abby 


Channel, 

Node 


Ridge deformation. 

Closed centers. 

Storm-storm (Ben, Carmen) 
Resultant vector 


Wavne 


Node, 

Channel 




Ike 


Channel, 
Closed centers 


Node 
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The node pattern was found in 1 1 of the 12 tropical cyclones. Only Typhoon 
Thad did not show the node pattern. Thad was unusual in that it displayed a broken 
ridge from the start. The break gave an excellent indication of the 72 h position in Fig. 
3.18. Official forecast errors were large due to numerical model predictions for 
strengthening of the ridge, w’hich did not verify (Joint Typhoon Warning Center, 1984). 
The channel pattern (Fig. 3.2a) was the next most frequently occurring pattern. 
Channels w'ere dominant in portions of the tracks of seven storms and secondary 
patterns in two other tropical cyclones. They were generally found in the low level field 
(700 mb) or in the late stages of a storm. The ridge deformation pattern (node 
intrusion, pinch, or break) w’as a dominant pattern in three storms (including Thad) 
and as a secondary pattern in four storms. The closed center pattern occurred almost 
as often as the ridge deformation pattern, but w’as a dominant pattern only in Typhoon 
Ike. The track of Ike (Fig. 3.19) is within the channel south of the zero contour tow’ard 
the closed +2 center before Ike turned toward the node to the west of Luzon. 

2. Cyclone-Cyclone Interaction 

The interaction between tropical cyclones is a complex problem for the 
forecaster. A determination first has to be made w'hether or not an interaction is 
occurring. The key questions are at what distance betw'een centers will the interaction 
begin and when will this interaction be strong enough to change either the track or 
intensity of one or both of the cyclones. According to Sandgathe (1987), storm- storm 
interaction takes place in four of the systems in the data set. STY Vanessa interacted 
to some degree with tw'o other tropical cyclones, but the track of Vanessa (not shown) 
does not appear to have been influenced during these interactions. The first interaction 
was a reduction of upper-level outflow above the incipient Vanessa due to deformation 
of the upper-level field by the outflow of TY Thad (Joint Typhoon Warning Center, 
1984). During the second interaction, Vanessa influenced the track of TY Warren, 
which was located in the South China Sea and had been tracking slowly westward. As 
STY Vanessa approached recurvature in the Philippine Sea, Warren turned rapidly to 
the east-northeast and increased in translation speed. This created a major forecast 
problem because of the potential danger to a group of ships leaving Subic Bay. 
Another storm-storm interaction between TY Clara and TY Bill resulted in little or no 
change in the track of Clara, but had a profound effect on Bill's track (not shown). 
The circulation and intensity of Clara at 12 UTC 18 Nov 1984 (Fig. 3.20) were very 
large compared to that of Bill. Clara's recurvature track was not changed by Bill, but 
Bill had a long track to the southeast under the influence of the circulation of Clara. 
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Fig. 3.18 Relative vorticity (xlO 's'^ in LAV at 12 UTC IS Oct 1984 
for Typhoon Thad. Forecast tracks are dashed with stars. 
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Fig. 3.19 Relative vonicity (xlO^s’^) in LAV at 12 UTC 30 Aug 1984 

for Typhoon Ike. 
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Fig. 3.20 Streamline analysis at 925 mb at 12 UTC IS Nov 1984 
illustrating Typhoon Bill and Typhoon Clara interaction 
(Joint Typhoon Warning Center, 19S4). 



Typhoon Gordon interacted with TY Faye at 00 UTC 30 Aug 82 (Fig. 3.21). 
The interaction is unusual because of the nearly 14 deg of latitude (840 n mi) distance 
at which it occurs. This distance seems particularly large considering that Faye is a 
very small cyclone. The interaction may be enhanced by the large positive vorticity 
(low) center northwest of Japan. 

In the fourth case of cyclone-cyclone interaction. STY Abby interacts with 
two other cyclones at e.xtended distances, and apparently before the other systems are 
even in warning status. Usually when cyclone-cyclone interaction is considered, two 
mature or nearly mature systems that are approximately 600 n mi apart are involved. 
ST^' Abby displays a very different pattern in the relative vorticity fields. About 36 h 
before the first warning on Tropical Storm (TS) Ben, an apparent interaction takes 
place between ST^' Abby (110 kt) and the tropical disturbance that would become TS 
Ben. The zero contour in the LAV field for 00 UTC 11 Aug 83 (Fig. 3.22a) surrounds 
Abby and a + 2 center to the east from which Ben will develop. The zero contour also 
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Fig. 3.21 Relative vorticiiy (xlO^s*^) in LAV at 00 UTC 30 Aug 1982 
illustrating Typhoon Gordon and Typhoon Faye interaction. 
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Fig. 3.22a Relative vorticity (xIO^s'^) in LAV at 00 UTC II Aug 1983 
illustrating Super Typhoon Abby - Typhoon Ben - Typhoon Carmen interaction. 
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Fig. 3.22b Relative vorticity (xlO^s’*) in LAV at 12 UTC 1 1 Aug 1983 
illustrating Super Typhoon Abby - Typhoon Ben - Typhoon Carmen interaction. 
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Fig. 3.22c Relative vorticity (xlO^s**) in LAV at 00 UTC 12 Aug 1983 
illustrating Super Typhoon Abby - Typhoon Ben - Typhoon Carmen interaction. 
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extends toward a + 2 center in the South China Sea that will spawn Tropical Storm 
Carmen 3 days hence. The pattern is further complicated by the strong -2 center to the 
north of Abby. In the next 12 h (Fig. 3.22b), the interaction between Abby and the 
disturbance in the South China Sea becomes stronger at a distance of approximately 15 
deg of latitude (900 n mi). The first available position for TS Ben is at 00 UTC 12 Aug 
(Fig. 3.22c) near 24.5°N, 145°E. In this field, interaction is taking place over a major 
portion of the western North Pacific basin and the South China Sea. The interaction 
may not just be multiple cyclone as large positive vorticity centers already exists that 
are adjacent (but separate from) both Ben and Carmen. The apparent interaction is 
taking place at a distance of 19 deg (1140 n mi) from the South China Sea to Abby 
and 16 deg (960 n mi) from Abby to the center to the northeast of Ben. Again, the 
strong ridge pattern north of STY Abby with the large negative (anticyclonic) center 
northeast of Ben further complicates the forecast situation. 

3. Forecast Performance 

Another consideration is which contours are giving the best overall forecast 
patterns and at what levels are the best forecast results found. Among the subjective 
aspects of the proposed technique are that different vorticity contours and different 
pressure levels may be involved as the storm evolves. The average rating values and 
time factors for the data set, by storm, are given in Table 8. Table 9 is a summary (by 
pressure level and by vorticity contour) of the occurrences, the average of the rating 
values and the average time factor (described above). Secondary node values are not 
given for the 700 mb level due to the small number of occurrences. Secondary nodes on 
the zero contour do not give high average ratings in this sample. The LAV field 
provides slightly better indications of track changes in this data set. This field may do 
even better if a mass w'eighted averaging scheme was used, as discussed previously. The 
-2 contour provides the most consistent information at all levels when it is present. 
Unfortunately, this contour is not as common as the zero contour in this data set. 

A test of the technique performance is how well it predicts the observed storm 
track turns. Table 10 is a compilation of the turn indication performance (hindcasts in 
this developmental sample) of the relative vorticity pattern technique. With the 
exception of the period of recurvature and northw'est movement for STY Abby, the 
technique appears very successful in this sample. Typhoon Gordon also had a vorticity 
pattern w'ith a stationary negative vorticity center to the east of Japan that is larger 
than the one that caused some problems in the forecast of STY Abby. In the case of 
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TABLE 9 

RELATIVE VORTICITY CONTOUR PERFOR.MANCE BY LEVEL 


Level 


Contour 


# of Fcsts 


Avq Rating 


Avg Time 


(mb) 

400 


2 


77 


(points) 
4.5 (good) 


(h) 

35 




0(p) 


169 


4.4 (fair) 


50 




0(s) 


78 


2.9 (poor) 


55 




-2 


87 


4.8 (good) 


58 


700 


2 


42 


4.7 (good) 


48 




0 


169 


4.3 (fair) 


49 




-2 


62 


4.5 (good) 


58 


LAV 


2 


60 


4.7 (good) 


38 




0{p) 


172 


4.6 (good) 


49 




0(s) 


67 


1 .9 (poor) 


44 




-2 


69 


5.1 (good) 


61 


(p) 


primary, (s) 


secondary 







TABLE 10 

TURN PREDICTION PERFORMANCE 
OF THE RELATIVE VORTICITY PATTERN TECHNIQUE 



Turn type 


Observed 


Forecast 


False Indicators 


1 Right 


8 


7 


0 


1 Left 


4 


4 


0 


1 Recurvature 


7 


7 


1 * 


Loop 




0 


0 



ST^' Abbv - Resultant-vector pattern sugcested recurvature while 
pnmar\' node indicated northwest movemen't. 



Gordon, recurv'ature would be forecast based on this technique. However, the pattern 
was not as clear as in the other storms, and the time factor was reduced. A clear 
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recurvature pattern did not appear until 24 h before the event. The experience with 
this set suggests that a blocking high is a definite complicating factor in forecasting 
recurvature. 

The one right turn that is not indicated by the synoptic pattern technique is a 
30 deg change in the track of TY Ike as the system turned toward the west under the 
subtropical ridge after moving to the southwest for approximately two days. This 
occurred at low latitudes (9°N) under conditions where recurvature was not a 
consideration. 

The loops in the the track of STY Abby and of TY Lex were short-term 
features (12 h or less), and were not evident in the 12 h time step of the fields. In any 
case, forecasting loops would probably cause the Typhoon Duty Officer's qualifications 
to be reevaluated and is not a high priority. 
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IV. CONCLUSIONS AND RECOMMENDATIONS 



This study indicates that relative vorticity charts have the potential to be a useful 
tool in the recurvature forecast decision process. Following the poor performance in 
the official forecasts of Super Typhoon Marge during October-November 1983, the 
JTWC initiated an internal review of the forecast process. This review did not reveal 
anything that could have been done with the existing maps and resources to improve 
the performance! This study of Marge suggests that a different presentation of the 
data available in the JTWC during the forecast development might have contributed to 
better forecasts. Thus, the first recommendation is that software be made available to 
JTWC to allow presentation of the Global Band Analysis fields in terms of relative 
vorticity contours. These vorticity fields allow the relative strengths and weaknesses of 
the major troughs and ridges to be analyzed readily, and incorporate dynamical 
features that are not apparent in the streamline presentation. A systematic procedure 
of locating the strong or weak points in the ridge may refine the recurvature longitude 
selection process. A mass-weighted LAV is suggested based on the experience with 
Marge. 

The developmental study had plotted best tracks on the relative vorticity contour 
charts to identify and locate the patterns described above. A true test of the 
operational importance and forecast capability of the technique should be carried out. 
The test should present an independent sample of storm sets in a manner similar to 
that found in the forecast center. The data presented should include only products and 
fields that would be available to the operational forecaster. Ideally, a first forecast 
should be developed for each 12 h period without having the relative vorticity fields 
available, followed at a later time by a second set of forecasts developed by the same 
forecaster having the relative vorticity fields. The skill level for each of the forecasts for 
a data set could be compared to a climatology-persistance forecast. The use of the 
same forecaster for each set would preclude individual inherent skills from biasing the 
results. To reduce any bias that would be introduced by the forecaster remembering 
tracks, a method of presentation that mixes the sequence of storms would be required. 
If the method helps eliminate from consideration forecast aids that are not likely to be 
correct, the translation speed may be better defined by the remaining aids. Further 
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study may also reveal a pattern in the relative vorticity fields that will help determine 
the translation speed. With the proper preparation and planning, the test could be 
carried out on a quasi real-time basis utilizing the data collection and display 
capabilities in the IDEA lab. 

The storm-storm interaction analysis as currently practiced in the forecast center 
could be refined through the use of the relative vorticity fields. The pattern seen in 
STY Abby and in the TY Gordon-TY Faye cases suggest that interaction may be 
taking place at greater than expected ranges and that complex interactions patterns 
may be involved. The ability to detect these features in the analyzed fields may give 
the forecaster a basis for evaluating guidance from a numerical model that considers 
only a single vortex in the field. The early detection of vortex interactions might alert 
the forecaster to these complex interactions. Thus, a study focused on the storm-storm 
interaction cases using the relative vorticity fields is recommended. 

The previous empirical orthogonal function (EOF) representations of storm- 
synoptic data by Wilson (1984), Schott (1985) and Weniger (1987) also might be 
applied to relative vorticity fields for eventual inclusion in the decision tree forecasting 
tool. A possible approach to stratification of the data set would be compositing of the 
data by the storm-synoptic parameters. The storm-synoptic approach could lead to 
improved forecast aid selection guidance. Another approach might be to composite 
the fields in a manner similar to the approach used by Allen (1984) in the COSVIOS 
technique. Application of relative vorticity information to COSMOS type field 
categories could lead to the development of an improved COSMOS technique. 

Future application and testing of this technique would be facilitated if the 
database in the IDEA lab included the best track data and the archived forecast aids. 
A study could then be undertaken to compare the relative vorticity patterns with the 
forecast aids for different categories of tropical cyclones. The objective is to select the 
"best" aid based on a set of rules developed from the past performance of the various 
aids under the conditions presented. The rules might be incorporated into an artificial 
intelligence shell on JTWC's computers for track aid selection. Even experienced 
forecasters could benefit from the detailed analysis of the new data required for 
inclusion in the rule base. Once the basic system is developed, the "experience" of the 
system could be continually updated as new information is available. This would 
provide a system for incorporating lessons learned in a particular forecast experience, 
which is now often lost when the forecaster leaves JTWC. Another advantage in the 
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development of an AI technique is that an interactive mode could be used in the 
training of new forecasters. Rules and decision points would be displayed to the trainee 
as a scenario develops. Appropriate and immediate feedback as decisions are made 
would develop the trainee's skills in a number of situations. 

Even though new dynamical models or other track guidance tools may be 
developed, the typhoon forecaster still must evaluate this guidance in relation to 
existing techniques. This study suggests that the relative vorticity fields would be a 
useful tool in displaying and understanding the interaction of the typhoon with 
adjacent synoptic features. It is expected that availability of these vorticity fields could 
assist the forecaster and hopefully improve track forecasts in the future. 
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APPENDIX A 

PROGRAM-GLOBAL BAND FILE IDENTIFICATION 

//SHGBID84 JOB (4697,9999)/SHERMAN 2200',CLASS = G 
//*MAIN LINES = 40 

//* 

/ i 

//* ** FILENAME FILETYPE ** 

//* GBFID NEDN 

I 

I i 

// EXEC FORTVCLG,PARM.FORT='LANGLVL(66)' 

//FORT.SYSIN DD * 

C 

C PGM TO IDENTIFY DATA FIELDS PRESENT IN GLOBAL BAND ANALS 
C ON FNOC TAPE 144X49 MERCATOR GRID NEDN FORMAT 
C 

C ORIGINAL AUTHOR UNKNOWN, UPDATED AND MODIFIED BY 
C B. T. SHERMAN, APRIL 1987 

C 

C THE PROGRAM AS PRESENTED WILL GENERATE A LIST OF THE DATA 
C FIELDS PRESENT ON A NEDN TAPE OF GBA FORMAT GRIDDED DATA. 

C TWO CHANGES ARE REQUIRED TO RUN THE PROGRAM. THE JOB NAME 
C AND USER ID INFO WILL NEED TO BE CHANGED TO GIVE A UNIQUE 
C IDENTIFIER TO THE OUTPUT. THE TAPE TO BE READ MUST BE 
C NAMED IN THE UNIT 10 SPECIFICATIONS AT THE END OF THE 
C PROGRAM 

REAL*4 ID(4745),BCD(23) 

INTEGER*4 IDTG(8),NP(8),ISKP(33),ITITLE(42),NC(3) 

INTEGER’^ INT(23),ITAU(4),BCDTAU(4) 

INTEGER IDATA(7056) 

INTEGER’='2 JDATA(14112) 

C 

DATA ENDS/'AAAAV 
C 

EQUIVALENCE (IDATA(1),JDATA(1)) 

C 

NA= 14220 
NR =7056 
NB = 0 
NCNT=0 
NSK = 0 
C 

CALL BUFFER(10,1,14220,ID) 

C 

C SKIP TO DESIRED DATA NSKPRC= NO. OF RECORDS SKIPPED 



95 



c 

READ (5,1) NSKPRC 

1 FORMAT (15) 

WRITE (6,2) NSKPRC 

2 FORMAT (' NO. OF RECORDS SKIPPED= ,15) 

IF (NSKPRC.EQ.O) GO TO 55 
IF (NSKPRCLE.IOOO) GO TO 5 
NSK = NSKPRC/ 1000 
D0 4 MSK=1,NSK 
CALL SKPREC(10,1000,&99,&999) 

4 CONTINUE 

NCNT= NSKPRC 

NSKPRC = NSKPRC-(NSK*1000) 

5 CALL SKPREC(10,NSKPRC,&99,&999) 

C 

55 CONTINUE 
C 

DO 926 L= 1,4745 
926 ID(L)=ENDS 
C 

CALL GETREC(10,&99,&9999) 

NCNT=NCNT + 1 
C 

CALL GBYTES(ID(1),NC,0,6,0,3) 

DO 31 J = 1,3 

31 INT(J) = NC(J) 

CALL GBYTES(ID(1),ITAU,24,6,0,4) 

CALL XBCD(4,ITAU,BCDTAU) 

CALL GBYTES(ID(2),IDTG,16,6,0,8) 

I = 1 

DO 32 J = 4,11 
INT(J) = IDTG(I) 

CALL XBCD(1,IN1(J),BCD(J)) 

32 I = I + 1 

C *** CONVERT BCD TO DECIMAL NUMBER •^***'^ 
CALL XINT(BCD,IYMDH,IER) 

CALL GBYTES(ID(4) ,NP,0.6,0,8) 

I = 1 

DO 33 J = 12,19 
INT(J) = NP(I) 

33 I = I + 1 

CALL GBYTES(ID(6),IUNIT,8,5,0,0) 

IN1(20) = lUNTT 

CALL GBYTES(ID(6),ISCLE,13,6,0,1) 

INT(21) = ISCLE 

CALL GBYTES(ID(7),M,24,12,0,0) 

IN 1(22) = M 
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CALL GBYTES(ID(8),X,4,12,0,0) 

INl(23) = N 

CALL GBYTES(ID(27),IDATA,8,16,0,NR) 

IF (lER.NE.O) GO TO 55 
CALL XBCD(20.IXLBCD) 

WRITE(6,203) IYMDH,(BCD(K),K= 1,20),NCNT 
203 FORMAT(1X/YMDH,BCD = MX,I10,2X,3A1,2X,8A1,2X,9A1, 

12X/RECORD NO. = ,18) 

GOTO 55 

9999 WRITE(06,120) 

120 FORMAT(lX; PARITY ERROR GETREC ') 

STOP 

99 WRITE(06,121) 

121 FORMAT(lX, ' EOF ENCOUNTERED IN GETREC ') 

STOP 

999 WRITE(06,122) 

122 FORMAT(lX/ SKPREC FAILURE ') 

STOP 

END 

&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& 

SUBROUTINE XINT(BCD, I YMDU,IER) 

REAL BCD(23),DIGITS( 10)/'0'/r/2','3'/4'/5'/6', 

#7'/8','97 

INTEGER INTGRS(10)/0,1,2,3,4,5,6,7,8,9/,INTS(8) 

DO 100 1 = 4,11 
DO 200 J= 1,10 

IF(BCD(I).NE.DIGITS(J))GO TO 200 
11 = 1-4+1 

INTS(II)=INTGRS(J) 

GO TO 100 
200 CONTINUE 

WRITE(6,1000) 

1000 FORMAT(2H ,'ERROR,CHECK YOUR SUBROUTINE XINT) 

IER=1 
RETURN 
100 CONTINUE 

IYMDH = INTS(l)* 10000000 

1 +INTS(2)* 1000000 

2 +INTS(3)*100000 

3 +INTS(4)* 10000 

4 +INTS(5)-1000 

5 + INTS(6)-‘T00 

6 +INTS(7)«10 

7 +INTS(8) 

IER = 0 
RETURN 
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END 

C &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& 

c 



c 

c 

c 

c 

c 

c 



SUBROUTINE XBCD(N, INPUT, BCDOUT) 



SUBROUTINE XBCD CONVERTS N INPUT CHAR.^CTERS FROM CDC 
EXTERNAL BCD TO IBM EBCDIC. THE INPUT CHARACTERS ARE 
STORED IN ARRAY INPUT AND THE OUTPUT (CONVERTED) 
CHARACTERS ARE STORED IN ARRAY BCDOUT. 



INTEGER HEX(47)/Z31,Z32,Z33,Z34,Z35,Z36,Z37,Z38,Z39, 
*Z21,Z22,Z23,Z24,Z25,Z26,Z27,Z28,Z29,Z12,Z13,Z14,Z15, 
Z16,Z17,Z18,Z19,Z0A,Z01,Z02,Z03,Z04,Z05,Z06,Z07,Z08, 

* Z09,Z30,Z20,Z10,Z11,Z3B,Z2B,Z2C,Z1B,Z1C,Z0B,Z0C/, 

* INPUT(N) 

r.EAL CHAR(47)/'A','B','C','D','E',T','G','H',T', 

* 'J','K','L','M','N','0',T','Q','R','S',T','U','V', 

* 'W','X','Y','Z','0',T','2','3','4','5','6',7','8', 

=«■ BCDOUT(N) 

DO 100 1 = 1, N 
DO 200 J= 1,47 

IF(INPUT(I).EQ.HEX(J)) GO TO 150 
200 CONTINUE 



J = 39 

150 BCDOUT(I) = CHAR(J) 
100 CONTINUE 



RETURN 

END 



r 

//LKED.SYSLIB DD 
// DD 
// DD 
// DD 
II DD 

II DD DISP = SHR,DSN = SYS3.TAPE.LOAD, LABEL = (,„IN) 

//GO.FTlOFOOl DD UNIT = 3400-5,VOL = SER= GB 1984, 

// DISP = (OLD,PASS), LABEL = (1,NL„IN), 

// DCB = (RECFM = F,BLKSIZE= 14220,DEN = 4) 

//GO.SYSIN DD 
0 



I* 

H 
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APPENDIX B 

PROGRAM-GET GLOBAL BAND ANALYSIS 

C THIS PROGRAM IS USED TO CREATE DATA TAPES FOR USE IN THE 
C IDEA LAB FROM NEDN FORMAT TAPES OF THE GLOBAL BAND ANAL. 

C THE PROGIU^M READS THE TAPE, CONVERTS THE NEDN FORMAT TO 
C IBM EBCDIC FORMAT FOR PROCESSING AND THEN WRITES THE 
C SELECTED FIELD DATA TO TAPE IN ASCII FORMAT. THE OUTPUT 
C FROM THE GBFID PROGRAM IS USED TO SELECT THE NUMBER OF 
C RECORDS TO BE SKIPPED TO ARRIVE AT THE START OF THE FIELDS 
C THAT COMPRISE A STORM DATA SET AND IDENTIFY THE RECORDS 
C THAT .MAKE UP THE DATA SET. IF RECORDS ARE MISSING FRO.M THE 
C GBA TAPES ENSURE THAT THE DATE/TIME INFOR FOR THE MISSING 
C FIELDS ARE NOT INCLUDED IN THE LIST OF DATA INPUT AT THE 
C END OF THE PROGRAM. SPECIFIC INFORMATION ON THE PAIU'\METERS 
C THAT WILL HAVE TO BE CHANGED FOR EACH DATA RUN IS SHOWN 
C AT THE END OF THE PROGRAM JCL STATEMENTS, NEAR THE LINES 
C TO BE CHANGED. REMOVE ALL CO.M.MENT LINES THAT ARE NOT IN 
C THE PROGRAM OR SUBROUTINES BEFORE ATTEMPTING TO RUN THIS 
C PROGRAM. 

C ORIGINAL AUTHOR UNKNOWN 

C UPDATED AND .MODIFIED BY B. T. SHERMAN, MAY 1987 

C 

//SG JOB ( ,9999),' ,2200',CLASS = G 

//*MAIN LINES = 40,RINGCHK = NO 

i Iff 
II 

ir 

// EXEC FORTVCLG,PAR.M.FORT = 'LANGLVL(66)' 

//FORT.SYSIN DD - 
C 

C PGM TO READ DATA FROM FNOC TAPE 144X49 MERCATOR GRID NEDN 
C FOR.MAT 
C 

DI .MENSION GLAT(49),G LON( 1 44) 

REAL-‘-4 ID(4745),BCD(23) 

REAL«4 DATA(7056) 

INTEGER*4 IDTG(8),NP(8),ISKP(33),ITITLE(42),NC{3) 

INTEGER-4 INl(23),ITAU(4),BCDTAU(4),TSTTWO/'2'/ 

INTEGER IDATA(7056) 

INTEGER*2 JDATA(14II2) 

C**NOTE** 

C DI.MENSIONS OF OUT MUST MATCH THE GRID SIZE OF OUTPUT 
C**NOTE** 

REAL OUT(I44,4I),D(I44,49) 
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c 

REAL XL(4) 

C 

DATA ENDS/'AAAA'/ 

DATA XL/'A'/E';G',T7 
CA-SFC E-700 G-400 T-250 1-200 
DATA NLVL/4/ 

DATA VRU/T7,VRV/'27 
C 

DATA C/.0436610743/,DD/ 114.59 15590262/ 
NAMELIST/NGRID/ NWST,NEST,NSTH,NNTH 
C 

EQUIVALENCE (IDATA(l)JDATA(l)) 

C 

NA = 28220 
NR =7056 
NB = 0 
NNS = 0 
NCNT=0 
INR = 0 
IOVER=0 
lER = 0 
A = 100. 

B = -100. 

CCC 

C 

C THESE ARE CURRENTLY FIXED, BUT CAN BE VARIABLE 
C 

CLAT = .1 
CLON = 179.9 

CCC 

C 

C DEFINE VALUES OF LAT,LON OF ALL GRID POINTS 
C 

C DO 2000 I =1,49 
C GLON(I) = 60. + (I-ir2.5 

C2000 CONTINUE 
C DO 2001 1 = 50,144 
C GLON(I) = ((I-49.)*2.5)-180. 

C2001 CONTINUE 
C DO 2002 11= 1,49 
C 1 = 50-11 

C GLAT(I) = 90.-ATAN(EXP(-ABS(II-3I.)*C))*DD 

C IF (II.GT.3I) GLAT(I) = -GLAT(I) 

C2002 CONTINUE 

C 

C 
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** READ IN GRID SIZE 

READ (5,NGRID) 

WRITE (6,NGRID) 

NT = NWST + NEST+l 
NJ = NSTH + NNTH+1 

IUCOMP = 20 
IVCOMP=20 
NRECU=0 
NRECV = 0 

READ IN NO. OF TAPES TO BE USED 

READ (5,78) NTAPE 

78 FORMAT (12) 

C DO 1000 ITP=1,7 
JUNIT=10 

CALL BUFFER(JUNIT, 1,14220, ID) 

** READ IN THE NO. OF TIME PERIODS TO BE INCLUDED IN THE 
STEP OF SKIPPING RECORDS AND SKIP TO DESIRED DATA 

READ (5,79) KTP 

79 FORMAT (13) 

READ (5,7) NSKPRC 
7 FORMAT (15) 

CALL SKIP (JUNIT,NSKPRC,IEND) 

MTP = 0 
IYMDH = 0 
224 MTP=MTP+1 



.-{t -it * 



READ (5,225,END= 123) IY,MTH,JD,JH,IREW 
225 FORMAT (I2,2X,3I2,1X,I1) 

IF (IREW.EQ.l) CALL REWIND (JUNIT) 
KYMDH = I Y^'H 000000+ MTH* 10000 + JD*100 + JH 
IF (lYMDH.EQ.KYMDH) GO TO 75 

IF (MTP.LE.KTP) GO TO 751 

READ (5,79) LTP 

KTP=KTP + LTP 

READ (5,7) NSKPRC 

NSKPRC = NSKPRC-2 

CALL SKIP (JUNIT,NSKPRC,IEND) 

IF(IEND.NE.O) GO TO 88 
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c 

751 CONTINUE 
927 DO 926 L= 1,4745 
926 ID(L) = ENDS 
C 

CALL GETREC(JUNIT,&88,&9999) 
GO TO 331 

88 IF (NTAPE.EQ.l) GO TO 99 
CALL REWIND(JUNIT) 

J UNIT = J UNIT +1 
CALL BUFFER (JUNIT, 1,14220, ID) 
READ (5,79) LTP 
KTP=KTP+LTP 
READ (5,7) NSKPRC 
CALL SKIP (JUNIT,NSKPRC,IEND) 
GO TO 751 
331 CONTINUE 



CALL GBYTES(ID(1),NC,0,6,0,3) 

DO 31 J = 1,3 

31 INl(J) = NC(J) 

CALL GBYTES(ID(1),ITAU,24,6,0,4) 

CALL XBCD(4,ITAU,BCDTAU) 

CHECK TO SEE IF IT IS A FCST FIELD IF SO IGNORE IT 
IF (BCDTAU(4) .EQ. TSTTWO) GO TO 751 
CALL GBYTES(ID(2),IDTG, 16,6,0,8) 

I = 1 

DO 32 J = 4,11 
INl(J) = IDTG(I) 

CALL XBCD(1,IN1(J),BCD(J)) 

32 I = I + 1 

C CONVERT BCD TO DECIMAL NUMBER 
CALL XINT(BCD,IYMDH,IER) 

IF(IER.NE.O) GO TO 751 
CALL GBYTES(ID(4) ,NP,0,6,0,8) 

I = 1 

DO 33 J = 12,19 
INl(J) = NP(I) 

33 1 = 1+1 

CALL GBYTES(ID(6),IUNIT,8,5,0,0) 

IN1(20) = lUNIT 

CALL GBYTES(ID(6),ISCLE,13,6,0,1) 

IN1(21) = ISCLE 

CALL GBYTES(ID(7),M, 24, 12,0,0) 

INT(22) = M 
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CALL GBYTES(ID(8),N,4,12,0,0) 

1X1(23) = N 

CALL GBYTES(ID(27),IDATA,8,16,0,NR) 

C 

C CDC HAS ONE'S COMPLEMENT, IBM HAS TWO'S COMPLEMENT FOR 
C NEGATIVE'S. 

DO 30 I = 1,NR 

30 IF (IDATA(I) .LT. 0) IDATA(I) = IDATA(I) + 1 

CHECK YYMMDDHH TO GET THE DATA WE NEED *”'=**'’ 

CALL XBCD(20,IXT,BCD) 

IF (lYMDH-KYMDH) 751,75,224 

*r* ^ 

75 CONTINUE 

IF(BCD(3).NE.VRV .AND. BCD(3).NE.VRU) GOTO 751 

SET UP TO UNPACK 16 BIT DATA POINTS. THUS JDATA WHICH IS 
-2 HAS ONE DATA POINT PER WORD. DATA IS IN THE RIGHTMOST 
16 BITS OF IDATA WHICH IS *4. THE SCHEME BELOW DOES NOT 
WORK IF THE NBTSDP = 20 OR 08. 



L = 0 

NR2 = NR'2 

IF (ISCLE .LT. 32) GO TO 40 
ISL= 2”(ISCLE-32) 

DO 17 I = 2,NR2,2 
L = L + 1 

IDUM = JDATA(I) 

17 DATA(L) = FLOAT(IDUM)/FLOAT(ISL) 

GO TO 45 

40 ISL = 2*-(32-ISCLE) 

DO 18 I = 2,NR2,2 
L=L+1 

IDUM=JDATA(I) 

1 8 D AT A( L) = FLOAT( I DUM)- FLO AT( I S L) 
45 CONTINUE 

TEST FOR LEVEL 



LVL = 0 

DO 500 K= 1,NLVL 
LVL=K 

500 IF(BCD(1).EQ.XL(K)) GOTO 505 

501 GOTO 751 
505 CONTINUE 

510 CONTINUE 
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511 CONTINUE 



C 

C DETERMINE GRID POINTS IS,IE,JS,JE FOR THE OUTPUT DATA 
C 

CALL FINDGR(GLAT,49,CLAT,I,NSTH,N'NTH,JS,JE) 

CALL FINDGR(GLON,I44,CLON',-I,NWST,NEST,IS,IE) 

IF (IS.LE.O) IS=IS+ 144 
IF(IE.GT.I44) IE= IE-144 
C 

C CREATE THE OUTPUT VECTOR OUT 
C IF JS.LE.O OR JE.GT.49, DATA AT THAT GRID POINT OF 
C VECTOR OUT WILL BE CONSIDERED MISSING 
C 

JD = JS-1 
DO 550 J=1,NJ 
IDD = IS-I 
IFLAG = 0 
JD = JD+I 

IF(JD.GTJE) GO TO 551 
IF (JD.GT.0.ANDJD.LE.49) GO TO 539 
IFLAG= I 

539 DO 540 1= I, NT 

IDD=IDD+ I 

IF (IDD.GT.I44) IDD= IDD-I44 
IF (IDD.GT.IE) GO TO 555 
IF (IFLAG.NE.O) GO TO 538 
OUT(I,J) = D(IDD,JD) 

IF (OUT(U).LT.B.OR.OUT(I,J).GT.A) lER = I 
GO TO 540 
538 OUT(LJ)= 99999. 
lOVER = I 

540 CONTINUE 

550 CONTINUE 

GO TO 570 

551 WRITE (6,552) JS,JE,JD 

552 FORMAT (' ERROR IN IDENTIFYING LAST LATITUDE LINE'/' 
1JS= ',I3,2X,'JE= ',I3,2X,'JD = ',13) 

STOP 

555 WRITE (6,556) IS,IE,IDD 

556 FORMAT (' ERROR IN IDENTIFYING LAST LONGITUDE LINE'/' 
HS = ,I3,2X,TE= ',I3,2X,'IDD= ',13) 

STOP 

570 CONTINUE 
C 

IF(BCD(3).NE.VRU)GO TO 800 
NRECU = NRECU+1 

WRITE(6,I0) NRECU,IYMDH,(BCD(K),K= 1,20) 
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I ,CLAT,CLON,IS,IE,JS,JE 

1 ,D(IS,JS),D(IE,JE),OUT(l,l),OUT(NI,NJ) 

2 ,IER,IOVER 

10 FORMAT (1X/NRECU = M4,I10,1X,3A1,1X,8A1 
l,lX,9Al,2X;CLAT=',F5.i; CLON = ',F6.1,2X, 

2TS = M3/ IE = M3/ JS = ',13/ JE = ',I3,4F7.l,I3,Il) 

CC WRITE (lUCOMP) NRECU,(BCD(K),K= l,3),IYMDH,OUT 
WRITE(IUCOMP,610)NRECU,(BCD(K),K=1,3),IYMDH 
610 FORMAT(I4,1X,(3(A1,1X)),I10) 

WRITE(IUCOMP,620)OUT 
620 FORMAT(10F8.2) 
lOVER = 0 
lER = 0 
GO TO 751 

800 IF(BCD(3).NE.VRV)GO TO 751 
NRECV=NRECV+1 

WRITE(6,11) NRECV,IYMDH,(BCD(K),K= 1,20) 

1 ,CLAT,CLON,IS,IE,JS,JE 

1 ,D(IS,JS),D(IE,JE),OUT(l,l),OUT(NI,NJ) 

2 ,IER,IOVER 

11 FORMAT (1X,'NRECV=', 14, 110,1X,3A1,1X,8A1 
1,1X,9A1,2X,'CLAT= ',F5.1,' CLON= ',F6.1,2X, 

2'IS= ',13,' IE= ',13,' JS= ',13,' JE= ',I3,4F7.1, 13,11) 

CC WRITE (IVCOMP) NRECV,(BCD(K),K= l,3),IYMDH,OUT 
WRITE(IVCOMP,630)NRECV,(BCD(K),K = 1,3),IYMDH 
630 FORMAT(I4,1X,(3(A1,1X)),I10) 

WRITE(IVCOMP,640)OUT 
640 FORMAT(10F8.2) 
lOVER = 0 
lER = 0 
GO TO 751 
CC STOP 
C 

775 WRITE (6,776) IYMDH,IX2,KCNT 

776 FORMAT('ERROR IN RECORD ARRANGEMENT'/'IYMDH = ',! 10,2X, 
l'IX2 = ',I10,'RECORD NO. = ',I5) 

STOP 

777 WRITE(6,150) lYMDH 

150 FORMAT(lX,'PGM ENDED IYMDH = ',I10) 

STOP 

9999 WRITE(06,120) 

120 FORMAT(lX,' PARITY ERROR GETREC ') 

STOP 

99 WRITE(06,121) 

121 FORMAT(lX, ' EOF ENCOUNTERED IN GETREC ') 

STOP 

999 WRITE(06,122) 
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122 FORMAT(lX/ SKPREC FAILURE ') 

123 STOP 

END 



A A A A A A A A 



>5: 5^. 



SUBROUTINE SKIP 



jJ'. ;J: 5l', il« ?:« * >;< 



^ jJ-: 3t? jJ; >Js 



❖ ❖ 



SUBROUTINE SKIP (JUNTT,NSKPRC,IEND) 

THIS SUBROUTINE SKIPS NSKPRC RECORDS FROM JUNTT 



IEND = 0 



WRITE(06,2) NSKPRC 

2 FORMAT ( O'/NO. OF RECORDS SKIPPED= M5) 
IF (NSKPRC.EQ.O) GO TO 6 
IF (NSKPRC.LE.IOOO) GO TO 5 
NSK = NSKPRC/ 1000 
D0 4 MSK=1,NSK 

CALL SKPREC(JUNTT,1000,&99,&999) 

4 CONTINUE 

NSKPRC = NSKPRC-(NSK’M 000) 

5 CALL SKPREC(JUNTT, NSKPRC, &99,&999) 

6 CONTINUE 

GO TO 123 
99 IEND= 1 
RETURN 

999 WRITE(06,I22) 

122 FORMAT(lX,' SKPREC FAILURE ') 

STOP 

123 RETURN 

END 



SUBROUTINE FINDGR ’*'’*'*’*'****=I'****'t'************- 



SUBROUTINE FINDGR (A,N,P,IFLAG,NG1,NG2,MS,ME) 

C 

C THIS SUBROUTINE DETERMINES THE GRID LINE CLOSEST TO 
C THE CYCLONE CENTER 

C AND RETURNS THE BOUNDARY GRID LINE NOS. 

C 

C DEFINITION OF PAIUXMETERS 

C A = ARRAY OF N VALUES OF LAT/LON 

C P= POSITION (LAT/LON) OF CYCLONE CENTER 

C NGl = NO. OF GRID POINTS ON WEST/SOUTH SIDE OF P 

C NG2 = NO. OF GRID POINTS ON EAST/NORTH SIDE OF P 

C MS = GRID LINE NO. ON WEST/SOUTH SIDE OF P 

C ME = GRID LINE NO. ON EAST/NORTH SIDE OF P 

C IN THE CASE THAT P IS EXACTLY IN BETWEEN TWO GRID LINES, 

C IFLAG DETERMINES WHICH GRID LINE WILL BE CONSIDERED THE 
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c 

c 

c 

c 



CLOSEST 

IFLAG= + 1 GRID LINE WITH HIGHER VALUE THAN P WILL BE 
CHOSEN 

IFLAG = -I GRID LINE WITH LOWER VALUE THAN P WILL BE 
CHOSEN 



DIMENSION A(N) 



IF(P.GE.-177.5) GO TO 25 
LI = 49 
L2 = 50 
RETURN 
25 DO 30 1= I,N 
II = I+I 

IF(ILGT.N) 11= I 

IF (P.LT.A(I).OR.P.GT.A(II)) GO TO 30 
DI = P-A{I) 

D2 = A(II)-P 
IF{DI-D2) 31,32,33 

31 MS=I-NG1 

ME=I+NG2 
GO TO 50 

32 IF (IFLAG.GT.O) GO TO 33 

GO TO 31 

33 MS=II-NG1 

ME=II + NG2 
GO TO 50 
30 CONTINUE 
WRITE (6,39) P 

39 FORMAT(TAILURE TO FIND LAT/LON CLOSEST TO CYCLONE 
ICENTER7' LAT/LON=',F8.2) 

STOP 

50 RETURN 



END 



SUBROUTINE XINT(BCD,IYMDH,IER) 

C 

REAL BCD{23),DIGITS(10)/'0','1','2','3','4','5','6', 
I'7','8','97 

INTEGER INTGRS(10)/0,1,2,3,4,5,6,7,8,9/,INTS(8) 
C 

IER = 0 
C 

DO 100 1 = 4,11 
DO 200 J= 1,10 

IF(BCD(I).NE.DIGITS(J))GO TO 200 
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11 = 1-4+ 1 

INTS{II)=INTGRS(J) 

GO TO 100 
200 CONTINUE 

WRITE(6,1000) 

1000 FORMAT(2H /ERROR, CHECK YOUR SUBROUTINE XINT) 

IER= 1 
RETURN 
100 CONTINUE 

IYMDH = INTS(1)*10000000 

1 +INTS( 2)- 1000000 

2 +INTS(3)-100000 

3 +INTS(4)* 10000 

4 +INTS(5)-1000 

5 +INTS(6)*100 

6 +INTS(7)-10 

7 + INTS(8) 

RETURN 

END 

&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& 



C 

C 

c 

c 

c 

c 



SUBROUTINE XBCD(N, INPUT, BCDOUT) 

SUBROUTINE XBCD CONVERTS N INPUT CHARACTERS FROM CDC 
EXTERNAL BCD TO IBM EBCDIC. THE INPUT CHARACTERS ARE 
STORED IN ARRAY INPUT AND THE OUTPUT (CONVERTED) 
CHARACTERS ARE STORED IN ARRAY BCDOUT. 



INTEGER HEX(47)/Z31,Z32,Z33,Z34,Z35,Z36,Z37,Z38,Z39, 
*Z21,Z23,Z24,Z25,Z26,Z27,Z28,Z29,Z12,Z13,Z14,Z15,Z16, 



* Z1 7,Z1 8,Z 1 9,Z0A,Z0 1 ,Z02,Z03,Z04,Z05,Z06,Z07,Z08,Z09, 

* Z30,Z20,Z10,Z1 1,Z3B,Z2B,Z2C,Z1B,Z1C,Z0B,Z0C/,INPUT(N) 
REAL CHAR(47)/'A','B','C','D','E',T','G','H',T','J' 

*'K','L','M'/N','0',T','Q','R','S',T','U','V','W', 

*'X','Y','Z','0','r,T/3','4','5','6','7','8','9', 

',7','.','S','=^-',')'/%','#'/@7,BCD0UT(N) 

DO 100 1= 1,N 



DO 200 J= 1,47 

IF(INPUT(I).EQ.HEX(J)) GO TO 150 
200 CONTINUE 



J = 39 



150 BCDOUT(I)=CHAR(J) 
100 CONTINUE 
RETURN 
END 



CCCC ALL OF THESE LIBRARIES ARE REQUIRED FOR THE PROGRAM TO 
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CCC WORK, THE PROGRAM ALSO CALLS A USER INACCESSIBLE ROUTINE 
CCC CALLED GETREC WHICH IS A PART OF THE TAPE UTILITIES 
CCC LIBRARY. ^REMOVE THESE FOUR LINES TO RUN* 

//LKED.SYSLIB DD 
// DD 

// dd 
// dd 
// dd 

// DD DISP = SHR, DSN = SYS3.TAPE.LOAD, LABEL = (,„IN) 

CCCC THESE NEXT TWO CONTROL STATEMENTS ARE THE READ AND 
CC WRITE FORMATS FOR THE TAPE UNITS. 

CC THE FIRST IS THE READ FORMAT 

CC FOLLOWING 'VOL = SER=' GIVE THE TAPE NAME(ON THE OUTSIDE) 
CC FILL IN DATA FROM THE TAPE SCAN AT THE APPROPRIATE SPOT 
CC LABEL = 

CC BLKSIZE = 

CC DEN = 

CC *****REMOVE THESE 9 COM.MENT LINES TO RUN *********** 
//GO.FTIOFOOl DD UNIT = 3400-5, VOL = SER = GB 1984, 



// 

II 

1 1 
1 1 



DISP = (OLD,KEEP),LABEL = (l,NL„IN), 

DCB = ( RECFM = F,BLKSIZE = 14220, DEN = 4) 

CC THE NEXT CONTROL STATEMENT IS THE WRITE TO TAPE STATEMENT 
CC THAT PUTS THE DATA ON THE TAPE IN ASCII FORMAT. 

CC FIRST CHANGE THE EXTERIOR LABEL LINE 'VOL = SER= XXXXXX 
CC NEXT YOU WILL NEED TO ADJUST THE BLKSIZE OF THE OUTPUT 
CC ONE FIELD SHOULD FILL ONE BLOCK FOR EFFICENT MANAGEMENT 
CC ON THE VAX, & BECAUSE THIS METHOD USES LESS TAPE/RECORD 
CC RECSIZE = (HEADER) + (8*GRIDSIZE) 

CC HEADER = 1*LRECL 

CC GRIDSIZE = NPTS(N7S)*NPTS(E/W) 

CC BLKSIZE MUST BE AN INTEGER MULTIPLE OF LRECL. 

CC IF RECSIZE IS NOT AN INTEGER MULTIPLE OF LRECL USE THE 
CC NEXT INTEGER TO GET BLKSIZE 
CC EX: IF RECSIZE = 63.4 

CC BLKSIZE = (64 * 80) = 5120 

CC ****REMOVE THESE 15 CO.MMENT LINES TO ruN***TCAD****** 
//GO.FT20F001 DD UNIT = 3400-5, VOL = SER= DSGB84, 



// DISP = (OLD,KEEP),LABEL = (l,NL„OUT), 



// DCB = (OPTCD = Q,RECFM = FB,LRECL = 80,BLKSIZE = 28320,DEN = 4) 

//* 

/ ; 



//GO.SYSIN DD * 

&NGRID NWST = 48, NEST = 95, NSTH = 10, NNTH = 30&END 
1 THIS IS THE NUMBER OF TAPES TO BE READ 



1 

3040 THIS THE NUMBER OF RECORDS TO SKIP 
84 6 912 0 YEAR,MONTH,HOUR, REWIND CODE:0= NO, 1 = YES 
84 610 0 0 



109 



APPENDIX C 

PROGRAM-GLOBAL BAND TAPE TO DISK 



program GBTTOD 
c 

c This program reads Global Band Analysis data from tape 
c and writes it to a file on disk for further processing, 
c The program uses QIO system services to read the data, 

c The subroutine (mtread) which actually calls the QIO is 

c located in disk$2: [metdept. tape]mt/lib. 
c (GBTTOD means Global Band, Tape TO Disk) 
c 

c Brett Sherman & Jim Cowie 7/87 (NFS) 

c 

c 

★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★ 

c 

c 

character data*28320 
integer*2 data_len, real_len 

integer*4 is tat, status, channel, sysSassign, sys$deassign 
parameter (data_len*28320) 
c begin 

c assign the QIO channel 

istat = sys$assign( 'tape', channel,,) 
write(*,*)' channel #', channel 
if(.not. istat)then 

write(*,*)' could””nt assign QIO channel number' 
write(*,*)' istat=' , istat 
goto 999 
end if 

c open the disk output data file 

open (unit=l, f ile*'gbanal.dat' ,status*'new' ) 
c using QIO's, read the record from tape and write to disk 

do 20 j=l,4 

100 call mtread(channel,%ref (data) ,data_len, status, real_len) 

write (*,*) ' real_len=' , real_len 
if ( .not. status) then 

call lib$signal(%val(status)) 
go to 999 
end if 

do 10 i=l, 28320, 80 
write(l,200)data(i: i+79) 
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10 continue 

20 continue 

c go back and read the next header and grid 
goto 100 

c deassign the QIO channel and exit the program 
999 is tat = sys$dassgn( channel ) 

200 format(a) 
stop 
end 
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APPENDIX D 

PROGRAM-FILE CONVERTER 



program FILE_CONVERTER 

cc This program converts the FNOC Global Band field data located in a 
cc disk file created by GBTTOD.FOR, to a gempak grid file, the program 
cc uses gempak subroutines to create the navigation data, the grid file, 
cc & write the data to the new file, each file contains a complete 
cc storm data set. 
cc 

cc Initial data relating to the grid size and lat/lon will have to be 
cc adjusted to values corresponging to the grid being converted. Ngrid 
cc will have to be changed to reflect the number of grids in the data 
cc set. The file names for read file, header file and the gem file will 
cc have to be changed for each different storm data set. 
cc 

cc In order for the grid conversion to run properly the LAT.DAT file 
cc must be accessible by the program. The program converts the 2.5 
cc degree longitude square GBA grid to a 2.5 deg latitude by 2.5 deg 
cc longitude grid for use in the GEMPAK routines, 
cc 

cc LCDR B.T. SHERMAN & JAMES COWIE JUNE 1987 

cc 

cc 

cc 

cc 



character file*80, rfile*80, setname*80, parm*4, lvl*l, hdrfil*80 
integer level(2) 

character time(2)*13, yymmdd*6, hh*2, mm*2, latfil*80 

logical replace 

integer fn, type, comp, grdno 

dimension grid(80,40) , rnvblk(256) , ighdr(2) , tgrid(7056) ,anblk(128) 
real intlat,igrid(144,49) ,mgrid(144,40) ,lat(49) 



c define the file to be read & the name of the file created 
c 

data rfile /'THAD_1984.DAT' / 
data file /'THAD 1984. GEM'/ 
data hdrfil /'THAD"1984.HDR' / 
data latfil /'LAT.DAT'/ 



c 

open (unit=4, file=latf il,status='old' ) 
open (unit=5, f ile=rfile,status='old' ) 
open (unit*6, file=hdrfil,status='new' ) 
c 

replace = .false. 
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c 

c read the latitude data 
c 

jl = 1 

do 4 i » 1,7 

read(4,5) (lat(j), j = jl , jl + 7) 

5 format(8(f9.3)) 

jl = jl -H 8 
4 continue 

c do 6 m = 1,49 
c write(6,*) lat(m) 
c 6 continue 

c 
c 

c set initial parameters 
c 

c ngrid is the number of grids in the data file to be converted, kx & 
cc ky are dimensions of the grid, the lat/log limit designators are set 
cc up as follows ex: rlatb = the real value of the latitude at the 
cc bottom of the grid, mm = the minutes for the time designator, ivcord 
cc = 1 designates pressure as the vertical coordinate, navsz, ihdrsz , 
cc ianlsz are set by the subroutines in the gempak, level(2) = -1 is the 
cc indicator that you are using only one level for the grid 
c 

ngrid = 104 
kx = 80 
igx = kx 
ky = 40 
igy = ky 
npts = kx * ky 
rlatb = -40.000 
rlatt = 57.500 
rlonr = 257.500 
rlonl = 60.000 
navsz = 256 
ianlsz = 128 
ihdrsz = 2 
ivcord = 1 
data mm /'OO'/ 
c 

c call the subroutine to create the nav data for the grid 
c 

call gr_mllnv (kx,ky, rlatb, rlonl, rlatt, rlonr, rnvblk,iret) 
c 

c if the return code indicates a problem print it and stop 
c 

if (iret.ne.O) then 

write(*,*) 'gr_mllnv error - iret=',iret 
go to 999 
end if 
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c call the subroutine to create the new data file 
c 

call gd_cret ( file , navsz , rnvblk, ianlsz , anlblk, file , ihdrsz , if Ino , i re t ) 
if (iret.ne.O) then 

write(*,*) 'gd_cret error - iret=',iret 
go to 999 
end if 
c 

grdno = 0 
c 

do 300 1 = l,ngrid 

read (5,10) fn,lvl, type , comp , yymmdd , hh 
10 format (IX, i3, lx,al, lx, il, lx, il ,3x,a6,a2) 

c 

WRITE ( * , 10 ) FN , LVL , TYPE , COMP , YYMMDD , HH 
grdno = grdno + 1 

wr i te( 6 , 15 ) f n , Ivl , type , comp , yymmdd , hh , grdno 

15 format(lx, i3, lx,al, lx, il, lx, il,3x,a6,a2, lx, i4) 
time(l) = yymmdd// W/hh//mm 

time(2) = ' ' 
c 
c 

if (Ivl.eq.'A') then 

level(l)=1000 
else if (Ivl.eq.'E') then 
level(l)=700 

else if (Ivl.eq.'G') then 
level( 1)=400 

else if (Ivl.eq.'T') then 
level(l)=250 
c 

end if 
c 

level(2)*-l 

c 

if (comp.eq.l) then 

parm »'UVND' 

else if (comp.eq.2) then 
parm »'VWND' 
end if 
c 

read (5,16) 

16 format(SOx) 
il=l 

do 400 m=l,706 

read (5,20) (tgrid(I), i =il,il+9) 

20 format (10(f8.2)) 

il=il+10 
400 continue 
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c 



vrite(6,*)' tgrid (1) »',tgrid(l) 



c 

c 

c 



kount = 0 
do 100 j = 1,49 

do 200 k « 1,144 

kount » kount + 1 

igrid (k,50-j) = tgrid (kount) 

200 continue 

100 continue 

interpolate the GBA grid to NOGAPS format 

do 510 i-1,144 

ragrid(i,17)»igrid(i,19) 

intlat*2.5 

intj»18 

: North of the Equator 

do 520 j=20,48 

if (lat(j ).lt.intlat.and.lat(j+l).gt.intlat)then 
ragrid(i,intj)=igrid(i, j)+(igrid(i, j+l)-igrid(i,j))* 
* (intlat-lat(j ))/(lat(j+l)-lat(j )) 

% 

intj » intj + 1 
intlat « intlat + 2.5 
end if 

520 continue 



c 

c 

c 



South of the Equator 

intj » 16 
intlat = -2.5 



do 530 j = 18, 2,-1 

if (lat(j ).gt. intlat. and. lat(j-l). It. intlat)then 
ragrid(i,intj)»igrid(i, j)+(igrid(i, j-l)-igrid(i, j))* 
* (intlat-lat(j))/(lat(j-l)-lat(j)) 

intj » intj - 1 
intlat = intlat - 2.5 
end if 

530 .. continue 

510 continue 



c make the grid to go to the GEM file 
c 

do 540 j = 1,40 

do 550 k = 1,80 

grid(k,j) » mgrid(k,j) 

550 continue 

540 continue 
c 
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c 

write(*,*)grdno, time, level, parm 



call gd_vrit ( if Ino, grid, igx,igy,ighdr, time, level, ivcord, parm, 
# nbsiz, replace, iret) 

if (iret.ne.O) then 

write(’^,*) 'gd_vrit error - iret=',iret 
go to 999 

end if 

300 continue 

999 stop 
end 



This latitude information is required by the grid conversion 
subroutine in the program FILE_C0NVERT (remove these three lines prior 
to processing data) 



-40.956 

-24.226 

-04.990 

14.821 

32.986 

48.069 

59.745 



-39.042 

-21.928 

-02.497 

17.222 

35.056 

49.711 

0.0 



-37.076 

-19.593 

00.000 

19.593 

37.076 

51.300 

0.0 



-35.056 

-17.222 

02.497 

21.928 

39.042 

52.835 

0.0 



-32.986 

-14.821 

04.990 

24.226 

40.956 

54.318 

0.0 



-30.866 

-12.393 

07.473 

26.484 

42.816 

55.750 

0.0 



-28.698 

-09.943 

09.943 

28.698 

44.621- 

57.131 

0.0 



-26.484 

-07.473 

12.393 

30.866 

46.372 

58.462 

0.0 
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